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Preface 


This report is the AFCL Users Guide to operating LOW!RAN 7. A more detailed scientific report 
describing the basic physical assumptions and theories utilized by the LOWTRAN models as well as the 
documented verification of the code will be forthcoming in a separate report. 

Section 1.1 contains a brief discussion of the LOWTRAN model followed by a description of the 
basic capabilities and the additional features added to the LOWTRAN 7 computer code, Scetion 2 
discusses the code structure, portability and availability. The complete operating instructions for 
LOWTRAN 7 are described in Section 3. Section 3 includes detailed explanations for the various nodes 
of execution of the program, the many choices of utilizing models (molecular, aerosols and or 
particulates) included in the code, as well as the different methods of replacing these models by user 
inserted data. 

A discussion of 4 test cases and a description of the program ovtput to tapes 6, 7, and 8 are 
included in Section 4. 

The three Appendixes describe the LOWTRAN 7 plot program, filter function program and newly 


written. scanning function program. 
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Users Guide to LOWTRAN 7 


1, INTRODUCTION 


1.1 The LOWTRAN 7 Model 


The LOWTRAN 7 model! and computer code ealeulates atinospherie transmittance and 
background radiance for a given atmospheric path at low spectral resolution, This version is an 
extension and update of the current code, LOWIRAN G2 (and its predecessors LOWTRAN 5°, 

LOWTRAN 44, LOWERAN 3B9, LOWTRAN 38, and LOWTRAN 27). All the options and capabilities of the 
LOWTRAN 6 code have been retained. 


(Received for Publication 15 Angust 1988) 
1. Kneizys, F X., Shettle, E.P., Abreu, L.W., Anderson, G.P., Chetwynd, J.tt.. Galery, W.O., Selby, 
JE.LA., aid Clough, S.A. (1989) Atrnosphene Transmittance / Radiance: The LOWTRAN 7 
Model, (in preparation). 
Knelzys, FLX., Shettle, E.P., Gallery, W.O., Chetwyna, JUL. Abreu, T.W., Selby, J.E.A., Clough, S.A, 
and Fenn, R.W. (1983) Atmospheric Transmillance/Radiance: Computer Code LOWTRAN 6, 
AFPGL-TR- 83-0187, (NTIS AD A137796), 
Kneizys, FX., Shettie, E.P.. Gallery, W.O., Chetwynd, J.H., Abreu, L.W., Selby, J.B.A., Pent, R.W,, 
and MeClatchey, R.A. (1980) AQnospheric Transmittance/Radiance: Computer Code 
LOWTRAN 5, AFGL-TR-80-0067, AD A0S8613. 
4A. Selby, J-B.A., Kneizys, F.X., Chetwynd, Jr., J... and MeClatchey, R.A. (1978) Atmospheric 
Transmittance/Radiance: Computer Code LOWTRAN 4, AFGI-TR-78-0053, AD A058643 

5, Selby, J.E.A., Shettle, E.P., and MeClatchey, R.A. (1976) Atmospheric Transmittance from O25 
10 28.5 um Supplement LOWTRAN 3B, AFGL-TR-76 0258, AD A0dO070 1, 

6G. Selby, LEA, and MeClatehey, R.A) (1975) Atmospheric Transmittance from 0.28 to 28.5 um 
Computer Code LOWTRAN 3, AFCISL-TR-75-0255, AD AQ17734. 


to 


Mw 


7. Selby, JLE.A. and McClatehey, R.A. (1972) Atmospheric Transmittance from 0.25 fo 28.5 wie 
Computer Code LOWTRAN 2, AFCRL-TR-72-0745, AD 763721. 


The LOWIRAN 7 code calculates atmospheric transmittance, agnospheric background tadiance, 
single scattered Solar and lunar radiuiee, direct solar drradianee, and unailtiple scattered solar and 
thennal radianee. The spectral resolution of the model is 20 cm ! (tall width ai hall maximum) in 
steps of 5 cn? from O to 50,000 em! (0.2 yum to infinity). A stile parameter band model is used for 
molecular line absorption and the effects of molecular continuumetype absorption; molecular 
scattering, aerosol, and hydrometeor absorption ang scattering are included. Retraction and earth 
curvature are considered in the caleulation of (he atmospheric slant path and attenuation amounts 
along the path. Representative atinospheric, aerosol, cloud, and rain models are provided in the code 
with options to replace them with user-provided theoretical or measured values. 

Anew atmospheric data base® consisting of separate molecular profiles (O to 100 kim) for thirteen 
(13) miuinor and trace pases is provided tor use with the LOWTRAN 7 model. Six referenee atmospheres, 
each defining temperature, pressure, density and mixing ratios for WoO, O4, CHL. CO, and NoO, allas a . 
funetion of alitade (selected from the U.S. Standard Supplements, 1966" and the U.S. Standard 
Atmosphere 197619) allow a range of climatological choices. - 

Separate band models and band model absorption parameters developed by Picrluissi and 
Maragoudakis!! are inchided in the LOWPRAN 7 model and code for the following molecules: HyO, Ox, 
NjO, CHy, CO, Oy, COg, NO, NO», NH, and SOg. Analytic causimission functions (double-exponential) 
replace numerical tables used in previous LOWPRAN models “hen aw models were developed with and 
based on degraded line-by-line spectral? and validated agai cet Jaboratory measurcimenis, 

The water vapor continuum absorption af 10 pum was msediied from that of LOWPRAN G. ‘The 
sell-density dependent continuum absorption values in this region were reduced approximately 20 
percent based on Bureu's labors tory measurements! 344 and the atmospheric measurcments of Devir 
et ay,!5 


8. Anderson, G.P., Clough, S.A., Kneizys, F.X., Chetwynd, J.EL, and Shietthe, E.P. (1986) APGL 
Atmospheric Constituent Profiles (0 120 kr, AFGL-TR-86-0110, AD A175173. 

9. NASA (1966) U.S. Standard Aunosphere Supplements, 1966, U.S. Government Printing Office, 
Washingion, DC, : 

10. NASA (1976) U.S. Standard Atmosphere Supplements, 1976, U.S. Government Printing Olfice. 
Washington, DC, 

1). Peirluissi, JW. and Maragoudakis, C.le. (1986) Molecidar Transmission Band Models for 
LOWTRAN, AFGL-TR-86-0272, AD A1L80655, 

12. Clough, S.A., Kneizys, FLX., Shetdle, .P., and Anderson, G.P. (1986) Attnospheric radiance and 
transmittance: FASCOD2, Proc. of the Sixth Conference on Atmospheric Radiation, 
Williamsburg, VA, American Meteorological Society, Boston, MA, 111-141, 

13, Burch, DE. and Alt, RL. (1984) Contin Absorption by E30 in the 700 1200 em ? and 2-400: 
2800 cnv? Windows, AFGL-TR-84-0128, AD AL47391. 

1.0 Bureh, DAE, (1985) Absorption by HO in Narrow Windows benecen 8000 and 4200 cnet, APGL- 

TR-85-00386, AD A166648. 


I>. Devir, A.D.. Ben-Shalom, A., Lipson, $.G., Oppenheim, ULPL. and Ribak. bE. (2985) Annospheric 
Transmittance Measurements: Compartson tottht LOWTRAN 6, Report RAA/99 85, Peehnmion- 
Israel Insitute cf Technology, Haifa 32060, Israel, 


New ultraviolet absurption parameters for molecular oxygen (Schumann-Runge bands, Herzbery 
continuum) have been added to the code, }G.17.18.19.20 The absorption data fer ezone in this region 
(Hartley and Huggins bands) has been updated or improved based om more recent data? !.22.23,24 
including the addition of Cemperature-dependent absorption coeflicients, 

An improved extra-terrestrial solar source function is tnueluded in the LOWLRAN 7 model. The 
data for this function are based on the work of Vantoosier et al? >26) Neckel and Labs27, Werhli28, 
and Thekeakara.2" It covers the speetral region from O to 57,470 em | and is generally compatible 
with the resolution of the molectdar absorption parameters of the LOWTRAN model. 


16. Yoshino, K., Cheung, A.S-C. Esmond, &.R., Parkinson, W.EL, Freeman, D.1E., Guberman, S.1L.,, 
Jenouvrier, A., Coquart, B., and Merienne, M.F. (1988) Improved Absorption Cross Sections of 
Oxygen in the Wavelength Region 205-240 mn of the Herzbery Continuum, Planet. Space Sci. 
36:1469-1475. 

17. Johnston, H.S., Paige, M., and Yao, F. (1984) Oxygen absorption cross sections in the Herzberg 
continuum and between 206 and 327K, J. Geophus, Res. 89:11661-11665. 


18. Shardanand (1977) Nitrogen: induced absorption of oxygen in the Herzberg continuum, J. Quant 
Spectrosc. Radiat. Transfer 18:525-530. 


19, Frederick, JE. and Hudson, RD. (1979) Predissociation Hnewidths and oscillator strengths for 
the (2-0) to (13-0) Schumann-Runge Bands of Oy, J. Molec. Spectrosc. 74:2.47-258, 


20. Yoshino, K., Freeman, D0... and Parkinson, Wil. (1984) Atlas of the Schumann Runge 
absorption bands of O» it: the wavelength region 175-205 um., J. Phys. Chem. 13:207, 


21.) Bass, A.M. and Paur, RuJ. (1985) The Ulraviolet Cross- Sections af Ozone. 1. Measuretaents, in 
Attmaspherie Ozone, Proceedings of the Quadrennial Ozone Symposium in Haucidiki, Greece, 
edited by C, Zeleros and A, Ghaz, pp. 606-616, D, Reidel, Inc, 

22. Molina, L.T. and Molina, M.d. (1986) Absolute Absorption Cross Sections of Ozone in the 
185-350 nm Wavelength Range, J. Geophys. Res, 91:14501-14509, 


23. Yoshine K., Freernan, D.E., Esmond, J.R., and Parkinson, W.H. (1988) Absolute Absorption Cross 
Section Measurements of Ozone in the Wavelength Region 238-335 ni and the Temperature 
Dependence, Planet, Space Sci. 36:395-398. 

ft. Caccianl, M., disarra, A., and Fioces, G. (1987) Laboratory Measurements of the Ozone 
Absorption Coefficients in the Wavelength Region 339- 362 nm at Different Te mperatures, 
Dept. of Physics, University of Roma - La Sapienza, llaly, Internal Nate No. 882, 


25. VanHoosier, M.E., Bartoe, J.D.. Brueckner, G.E., and Prinz, D.K. (1988) Absolute Solar Spectral 
Irradiance 120 nm-400 nim (Results from the Solar Ultraviolet Spectral Irradiance Monitor- 
SUSIM-Experiment on Board Spacelab 2, Astro, Lett, and Commuiications 27: 163-168. 

26. Vantoosier, MLE. and Broecker, GE. (1987) Solar Ultraviolet Spectral Irradianee Monitor 
(SUSIM): Calibration Results from Spacel.b 2. Proceedings of the 8h Workshop on Vacuum 
Uuraviolet Radiometric Calibration ef Space Experiments, 19 March 1987. 

27. Neekel, Hand Labs, D. (1984) The solar radiation between 3300 and 12500 A, Solar Phys, 
90:2.05-2.58, 

28. Webhrli, Ch. (1985) Eeatra terrestial Solar Spectrum, Publication No, 615, July 1985, 
Physikalisch-Meteorologisches Observatorium and World Radiation Center, CH-7260 Davos 
Dorf, Swilzerland. 


29. Thekeakera, M.P, (1974) Extraterrestrial Solar Spectrum, 83000-6100 A at 1 A Intervals, Appl. 
Opt. 13:518-522. 


An efficient and accurate mulliple scattering parameterization3?.31 has been implemented in 
the LOWTRAN model based on the two stream approximation and an adding method for combining 
atmospheric layers, An interface scheme was also developed using the k-distribution method to 
match the multiple scattering approach to the LOWTRAN band model calculation of molecular gaseous 
absorption. The error of the multiple scattering parameterization in solar and themnal radiance 
calculations considering all possible viewing angles is estimated to be less than 20 percent. 

ror LOWTRAN 7, all the existing aerosol models and the rain model in the previous LOWTRAN 6 
model were extended through the millimeter wavelength region. In addition, the Navy Maritime 
model was modified to improve its wind speed dependence for the large particle component.* Water 
cloud models (cumiulus, stratus, altostratus, strato-stratocumulus and nimbostratus) from 
FASCOD2!2 have also been added. Two new cirrus cloud models with a realistic wavelength 
dependence and separate absorption, scattering and asymmetry parameters were developed for 
LOWTRAN 7 as well as a new aerosol model for desert conditions with a wind speed dependence.32:33 
The prograrn now provides for modified aerosol profiles over elevated surfaces, commonly referred to 
as the "Denver" case.34 

For the stratospheric aerosols, additional combinations of the wavelength dependent extinction 
coefficient models (background stratospheric, aged volcanic or fresh volcanic) and the vertical 
distribution profiles (background and moderaie, high or extreme volcanic), are available. The 
background stratospheric extinction model has been modified to utilize new refractive index data aud 
sive distribution measurements.3° 


2. DISCUSSION OF LOWTRAN 7 CODE 


This section describes the LOWTRAN 7 program structure, Six tables are included that contain 
names aud descriptions of each subroutine and block data 1sutine as they appear within the program 
structure. Section 2.1 presents five figures depicting the program structure. Section 2.2 discusses 
program portability, FORTRAN compatibility, precision, and some specific comments related to 


20. Isaacs, R.G., Wang, W-C., Worsham. R.D., and Goldberg, S. (1986) Multiple Scattering Treutmeni 
for use in the LOWTRAN and FASCODE Models, AFGL-TR-86-0073, AD A173990. 


31, Isaacs, R.G., Wang, W-C., Worsham, R.D., and Goldenberg, S. (1987) Multiple Scattering 
LOWTRAN and FASCODE Models, Appl. Opt. 26:1272-1281, 

32. Longtin, D.R., Shettle, E.P., Hummel, J.R., and Pryce, J.D. (1988) A Desert Aerosol Model for 
Radialive Transfer Studies, pg. 261-269 in Aerosols and Climate, ed. by P.V. Hobbs and M.P. 
McConmick, A, Deepak Publishing, Hampton, VA. 


33. Lonytin, D.R., Shettle, E.P., Hummel, J.R., and Pryce, J.D. (1988) A Wind Dependent Desert 
Aerosol Model: Radiative Properties, AFGL-TR-38-0112. 


34. Shettle, E.D. (1989) Comments on ihe Use of LOWTRAN in Transniission Calculations for Sites 
with Ground Level not at Sea Level, to appear in Appl. Opt. 28, 


35. Hummel, J.R., Shettle, E.P., and Longiin, D-R. (1988) A New Backgrouncl Stratospheric Aerosol 
Model for use in Atmospheric Radiation Models, AFGL-TR-88-0166. 


* Gathman, 5S, (1988) Private Communication. 


implementation of the program. Section 2.3 contains pertinent information on the program mn 
execution field length, and lastly Section 2.4 describes the availability of the program package. 


2.1 LOWTRAN 7 Code Structure 


A graphical representation of the LOWTRAN 7 main program structure is depicted in Figure 1. 
The structure for non-standard model, air mass, single scatlering geometry, and transmittance 
modules are shown in Figures 2 through 5 respectively. Descriptions of the major executable 
subroutines shown in Figure 1 are given in Table 1. The subroutines in the non-standard model, air 
mass, single scattering geometry and transmittance modules are explained in Tables 2 through 5 
respectively. Table 6 contains a brief description of the Block Data subroutines. 


2.2 Portability 


LOWTRAN 7 was developed on a CDC CYBER (a 60 bit-per-word machine) with the FORTRAN 77 
(FORTRAN 5) compiler. A major effort has been made to give LOWTRAN 7 the capability to run in 
singve precision on a 32 bit-per-word computer. 

LOWTRAN 7 uses four files: 

1, INPUT, read on UNIT = 5, containing LOWTRAN input directions. Maximum record length is 
80 characters. 

2. OUTPUT, written on UNIT = 6, containing the standard LOWTRAN output. 

3. TAPE7, written on UNIT = 7, containing copies of the input cards and the spectral results 
(transmittance and/or radiance). Used as input by plot, filter and scanning function programs (sce 
app2ndixes). 

4, TAPES, written on UNIT = 8, containing the transmittance for each individual gas. When 
running in the radiance mode the user may direct additional output to TAPE8 in one of the folk wing 
ways: 

a. The black body and differential transmittance layer values for the line of sight 
or when running with multiple scattering 
b. | The thermal and soiar vertical fluxes by layer. 


The unit numbers for these files are stored in the variables IRD(= 5), IPR(= 6), IPU(= 7), IPRI{= 8), 
which are carried in the common block IFIL. These files are accessed using OPEN slalements in the 
main program, 

bor computers with virtual nemory, the program must be compiled with "GLOBAL SAVE" or the 
local equivalent. Otherwise, information relating to the phase function is initialized in a first call to 
a subroutine, but lost in subseqent calls. 

Vie cirrus cloud rnodel includes an option to generate cirrus clouds at random altitudes. This 
option calls the machine dependent random tumber generator function subroutine RANF, which is 
called from RANDOM. The user will have to supply the local equivalent of RANF. 

The user may find it useful to have the date and {ime printed at the beginning of each case, The 
statements required to do this are commented out in the main program with 'C@' in columns 1 and 2. 
The subroulines FDATE and FCLOCK, which return the date and iine in A8 format, are commented 
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Table 1. Description of LOWTRAN 7 Subroutines 


LWTRN 7 - Main driver program, Reads control cards. 


MARINE -  Detennines aerosol extinction and absorption coefficients for the Navy 
maritime model. 
rn RDNSM - Reads user input data when MODEL 7 and VSA options are selected. 
JOU - Interpretive routine for JCHAR. 
. CHECK - Units conversion for pressure and temperature. 
DEFALT - Chooses a stored atmospheric profile and interpolates default values for 
a specific altitude. 
CONVRT - Accommodates uniform data input for mudel 0 o: 7. 
WATVAP - Computes water vapor number density (mol cm-) to accommodate 


"JCHAR" definitions specified by user. 


RDEXA - IHAZE7 or ICLD 11 triggers up to 4 regions of user input. 
AERNSM - Defines model atmosphere, aerosol profile and cloud profile. 
CIRRUS - Generates altitude profiles of cirrus cloud density. 


| RANDOM - Calis machine-dependent function RANF, that is a uniform random 
number generator. 
VSA - Army vertical structure algorithm of aerosol extinction and relative 
humidity for low visibility/low ceiling conditions. 
SUBSOL -  Caiculates the subsolar point angles based upon time and day. 
STDMDL - Sets up scaled densities from the model atmosphere. 
GEO - Driver for air mass subroutines. Calculates attenuator amourits for the 
slant path. (GEO can be called by both SSGEO and the main driver.) 
5 AERTMP - Defines temperature for each aerosol altitude region. 
SSGEO - Obtains attenuator amounts from scattering points along optical path to 
, the extraterrestial source. 
EXABIN - Loads aerosol extinction, absorption and scattering. 
EWULWC - Calculates liquid water conteni of standard aerosols, 
TRANS - Calculates transmittance, atmospheric radiance, and solar/lunar 


scattered radiance for slant path. Sets up data for double exponential 


band model. Evaluates vertical profiles of optical quantities required 


for mulliple scattering calculations. 


Table 2. Descripticn of AERNSM Subroutines 


ee 


AERNSM Defines model atmosphere, aerosol profile and cloud profile. 


CIRR18 - New cirnuis profile and new optically thin cirrus profile, 

FLAYZ - Final LOWTRAN altitude boundartes. 

JOU -— Interpretive routine for JCHAR. 

CHECK - Units conversion for Pressure and Temperature. 

DEFALT - Chooses a stored almospheric profile and intcrpolates default values 


for a specific altitude. 


CONVRT ~— Accommodates uniform data input for MODEL 0 or 7. 
WATVAP — Computes water vapor number density (mol cm) tu accommodaic 

“JCHAR" definitions specified by user. 

VSANSM - Used with VSA and MODEL 7, 

LAYVSA - Layering of atmosphere with VSA model. 

LAYCLD — Layering of atmosphere with ICLD 1 to 11. 

MARINE - Determines aerosol extinction and absorption coefficients for the Navy 
maritime model. 

DESATT - Determines aerosol extinction and absorption cvelficients, 

CLDPRF ~ Standard cloud profiles. 

AERPRE - Compute density profiles for aerosols. 


. | GEO 

: GEOINP 
FNDHMN 
REDUCE 


FDBETA 
RF PATH 


FILL 
LAYER 
| RADREF 


FINDSH 


SCALHT 
ANDEX 


‘ | | EXPINT 
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Table 3. Description of AIR MASS Subroutines 


Driver for air mass subroutines. Calculates attenuator amounts for 
the slant path 

Interprets geometry input parameters into the standard fonnH11, H2, 
ANGLE, and LEN. 

Calculates HMIN, the minimum altitude along the path and PHI, the 
zenith angle at H2. 

Eliminates slant path segments that extend beyond the highest profile 
altitude. 

Calculates angle, given H1, H2, and BETA by iicration. 

Determines the refracted path and the absorber amounts through all 
the layers. 

Defines the boundaries of the slant path and interpolates densities at 
these boundaries. 

Calculates the path and amounts through onc layer. 

Computes radius of curvature of the refracted ray for a horizontal path. 
Finds layer boundaries and scale height at ground for index of 
refraction. 

Calculates scale height of index of refraction. 

Computes index of refraction at a specific height. 

Perlorms exponential interpolations for the geomeiry routines, 


Teble 4. Description of SSGEO Subroutines ‘ 


SSCEO - Obtains attenuator amounts from seattering points along optical path 
to the extraterrestrial source. 
PSIDEL - Calculates the relative azimuth between the line of sight and the direct 
solar/lunar path. 


PSI -  Retums solar azimuth relative to line-of-sight at current scattering 
location, 

DEL - Returns solar zenith angle at any point along optical path. 

GEO - Driver for air mass subroutines. Calculates attenuator amounts for 


the slant path. 


SCTANG - Returns the scattering angle at any point along the optical path. 


TRANS 


MAPMS 
CXDTA 
C4DTA 
ABCDTA 
SLF296 
SLF260 
FRN296 
SINT 
C6DTA 


C8DTA 


HNO3 
AEREX7 


GAMFOG 


| Inpex 
| DEBYE 
HERIDA 


SCHRUN 
OSHHTO 


O3INT 
O3HHT] 
| O3INT 


Table S. Description of TRANS Subroutines 


Calculates transmittance, atmospheric radiance, and solar/lunar 
scattered radiance for siant path. Sets up data for double exponential 
band model, Evaluates vertical profiles of optical quantities required 
for multiple scattering calculation, 

Mapping routine from line of sight to vertical path. 

Locates coefficient for double exponential. 

Retums Ng continuum absorption coefficient at required wavenumber. 
Moves double exponential coefficients into new arrays. 

Loads self-broadened water vapor continuum at 296°K. 

Loads self-broadened water vapor continuum at 260°K. 

Loads foreign-broadened water vapor continuum at 296°K,. 

Performs interpolation for water vapor contint:um 

Returns Rayleigh molecular scattering attenuation coeflicient at required 
wavenumber. 

Returns Chappuis ozone visible absorption coefficient at required 
wavenumber (13000 to 24200 cn}). 

Determines nitric acid absorption coefficient at required wavenumber. 
Interpolates aerosol attenuation coefficients and asymmetry paranieter 
to required wavenumber. 

Computes attenuation of equivalent liquid water content in clouds. 
Calculates real and imaginary part of refractive index of water. 
Calculates wavenumber dependence of dielectric constant of water. 

UV Og Herzberg continuum - analytic funciion. 

UV Oy Schumann-Runge band model parameters, 

UV Og Hartley band temperature dependent coefficient: constant tenn 
(24370 ta 40800 enr}). 

Interpolation for Hartley band. 

UV Og artley band temperature dependent coefficient: Hnear term, 
Interpolation for Hartley band. 


Table 5. Description of TRANS Subroutines (Continued) 


OSHIT2 - UV Og Hartley band temperature dependent coefficient: quadratic 
tem. 
O3INT - Interpolation for HarQey band. * 
C3UV - UV Oz interpolation for 40800-54054 cnr ! region. 
O2CONT - Q continuum for 1395-1760 cnr? region. . 
| O2INT - Interpolation for Og continuum. 
TNRAIN - Calculates extinction due to rain as a function of rain rate. 
| GMRAIN - Computes attenuation of condensed water in the form of rain. 
RNSCAT - Calculates extinction, scattering and asymmetry p rameter due to 


rain in microwave region. 
DBLTX - Transmittance from new double exponential band model, 
SSRAD - Performs the layer by layer single scattering radiance sum, 
| PHASEF 


Chooses correct phase function based on relative humidity, frequency, 
scattering angle, and model. 
INTERP - Performs linear or logarithmic interpolation. 


PE -  Retums the appropriate phase function from the stored database. . 
HENGNS - Calculates phase function using Henyey-Greenstein method, 
MSRAD - Sets up profiles of optical properties for vertical path, Evaluates path 
integral of source function, multiple scattered radiance contribution, 
SOURCE - Contains solar intensity data and calculates lunar intensity. F 
FLAADD - Calculates upward and downward fluxes. Multiple scattered source a 
function is evaluated from fluxes using the stream approximation, 7 
SOURCE - Calculates solar intensity data and calculates lunar intensity. 
SUN - Evaluates extraterrestrial solar irradiances. 


MLATMB 
TITLE 
PRFDTA 
EXTDTA 
SF2Z96 
SF260 
BFH2O 
C4D 


MARDTA 
PHSDTA 
MDTA 
DSTDTA 


ATMCOM 
ABCD 


BOSHHO 
BOSHH1 
BOSHH2 
O3U0VFD 
BO2C 
CPRTRCG 
CPUMIX 
CPH20 
CPO3 
WVBNRG 
SHUMG 
SOLAR 


Table 6. Yescription of Block Data Subroutines 


Model atmospheric data. Six stored atmospheric models. 

Titles for output. 

Aerosol profile data. 

Acrosol and cloud extinction, absorption and asymmetry paranielers. 
Self-broadened absorption coefficients for water vaper centimuum at 296°K. 
Self-broadened absorption coefficients for water vapor continuum at 260°K, 


Foreign-broadened absorplion coetlicients for water vapor continuum at 296°K. 


Nitrogen continuum absorption coelficients and visibic vzone abscrption 
coefficients. 

Navy marine aerosol extinction and absorption data. 

70 averaged phase functions and truth table identifying correct phase fianetion. 
Cloud and rain modeled atmospheric data. 

Desert aerosol extinction, absorption coefficients and asynimetry parameters 
for 4 wind speeds. 

Inilializes constants used ii program. 

Stores the absorber model parameter for cach gas for the double exponential 
formulation and the coeflicients of each pas for the k-distribution. 
Contains Og Hartley-Huggins cross sections for 273K: constant term, 

Og Hart ley-Huggins cross sections: linear term. 

O3 Hartley-Huggins cross sections; quadratic terni, 

Og UV absorption coefficients (40800 to 54054 cn 1), 

Og continuum equivalent coefficients (1395-1760 cm=}). 

Band model absorption coefficients far trace gases, 

Band model absorption coefficients for uniformly mixed gases, 

Band model absorption coetucients for water vapor. 

Band model absorption coellicients for ozone. 

Specification for wave number band region limits for each gas absorber. 
Schinann - Runge Oy band model. 


Extraterrestrial solar spectral irradiances. 


—— 


and located at the ena of the main program LWTRN7. The user will have to modify these subroutines 
as necessary. 

The sample output included on the LOWTRAN 7 tape was generated on a CDC CYBER that has 14 
decimal digits of precision. The output for these same cases generated in a 32 bit-per-word machine 
with about seven decimal digits uf precision (for example, an IBM 370 or 4341 or a VAX) should agree 
with the sample to within about four decimal digits. When calculating the radiance where the 
absorber amounts are very small, the truncation error of a 32 bit-per-word machine might cause the 
radiance to be zero. 


2.3 Execution Field Length 


There are several techniques available to the user to reduce the executable field length of a 
program, On CDC CYBER sysiems, the program may be run using the SEGMENT loader, which 
effectively creates an overlay structure. File 2 of the LOWTRAN 7 tape (Section 2.4) gives the segment : 
loader directives that reduce the field length to less than 207,000,. Users of other systems who wish to 
create overlays should consult the structure chart in Figure 1 and the SEGMENT input directives for 
guidance. If certain program options are not required, the field length can be reduced by not including 
the corresponding subroutines. 


2.4 availability 


The LOWTRAN 7 package is available from: 


National Climatic Data Center, NOAA 

Environmental Data Services 

Federal Building 

Asheville, NC 28801 

(704) 259-0682 

The package is normally distributed on magnetic tape with the following characteristics: 

1. 9 track, 1600 BPI 
2. unlabcled 
3. ASCII 
4. fixed-length records, 140 characters-per-record, one record per block. 


The tape has 16 files. The contents of these files are as follows: 


LOWTRAN 7 source code 

SEGMENT input direction (relevant to CDC/NOS/BE systems only) 
INPUT file for test cases 

OUTPUT file for test cases 

TAPE 7 file for test cases 

TAPE & file for test cases 

Source code for plot program (see Appendix A) 

INPUT file for test cases of plot program 


CaAnNnNoa kwh 


OUTPUT file for test cases ef plot program 


° 


Source code for filter program (see Appendix B) 
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11. INPUT file for test cases of filter program 

12, OUTPUT file for test cases of filter program 

13, Source code for scanning function program (see Appendix C) 
14. INPUT file for test case of scanning function program 

15. OUTPUT file for test case of scanning function program 

16. List of all LOWTRAN phase functions. 


The record length for all the files is 140 characters; however, for files containing source code, 
INPUT files, and phase functions. only the first 80 characters contain infonnation, Users who 
require a different format should centact the National Climatic Data Center. The tape presently costs 
$101.00, 


3. INSTRU “TIONS FOR USING LOWTRAN 7 


The instructions for using LOWTRAN 7 are similar to those for previous versions. However, 
many new parameters have been added necessitating the addition of five new optional cards. The new 
parameters are principally required for reading in user supplied vertical! profiles of atmospheric 
constituent gases, rain and cloud inputs, and aerosols. 

In general, for standard atmospheric models, five input cards are required to run the program for 
a given problem. For a specific problem a combination of several of the fourteen additional optional 
control cards are possible. The formats for the five main cards, fourteen optional cards, and 
defisutions of the input parameters are given below. 


3.1 Input Data and Formats 


The use of the word 'CARD' is equivalent to editing with 80 columns. 
The program is activated by submission of a five (or more} card sequence as follows: 


CARD 1: MODEL, ITYPE, IEMSCT, IMULT, M1, M2, M3, M4, MS, M6, MDEF, IM, NOPRT, 
TBOUND, SALB 
FORMAT (1315, F8.3, F7.2) 


CARD 2; IHAZE, ISEASN, IVULCN, ICSTL, ICLD, IVSA, VIS, WSS, WH, RAINRT, 
GNDALT 
FORMAT (615, 5F10.3) 


OPTIONAL CARDS 


CARD 2A: CTHIK, CALT, CEXT, [SEED (If ICLD=18, 19, or 20) 
FORMAT (3F 10.3. 110) 

CARD 2B: ZCVSA, ZTVSA, ZINVSA (if IVSA=1) 
FORMAT (3F 10.3) 

CARD 2C: ML, IRD], IRD2, TITLE (if MODEL=0 or 7, and IM=1) 
FORMAT (315, 18A4) 


CARDS 2C1 THROUGH 2C3 (AS REQUIRED) ARE REPEATED ML TIMES 


CARD 2C1; ZMDL, P, T, WMOL(1), WMOL{2), WMOL(3), JCHAR 
FORMAT (F 10.3, 5E10.3, 15A1) 


CARD 2C2: (WMOL{J), J=4, 12) (If IRD 1=1) 
FORMAT (8E 10.3) 


CARD 2C3: AHAZE, EQLWCZ, RRATZ, IHA1, ICLD1, IVUL1, 


ISEAL, ICHR1 (If IRD2=1) 
FORMAT (10X, 3F 10.3, 515) 
CARD 2D: —IREG (1 TO 4) (If IHAZE=7 or iCLD=11) 
FORMAT (415) 


CARD 2D1: AWCCON, TITLE 
FORMAT (E10.3, 18A4) 


CARD 2D2: (VX(D, EXTC(N,D, ABSC(N,D, ASYM(N,]), l=1, 47) 
(if IHAZE=7 or ICLD=11) 
FORMAT (3(F6.2, 2F7.5, F6.4)) 


CARD 3: H1, H2, ANGLE, RANGE, BETA, RO, LEN 
FORMAT (6F 10.3,15) 


ALTERNATE CARD 3: H1, H2, ANGLE, IDAY, RO, ISOURC, ANGLEM (If IEMSCT=3) 
FORMAT (3F 10.3, I5, 5X, F10.3, 15, F10.3) 


OPTIONAL CARDS: 


CARD 3Al1: IPARM, IPH, IDAY, ISOURC (IEMSCT=2) 
FORMAT (415) 

CARD 3A2:; PARM1, PARM2, PARM3, PARM4, TIME, PSIPO, ANGLEM, G 
FORMAT (8F 10.3) (If IEMSCT=2) 


CARD 3B1: NANGLS (if IPH=1) 
FORMAT (15) 
CARD 3B2(1 TO NANGLS): (if IPH=1) 


(ANGF (), FQ), F(2.0, F(3,0, F(4., l=1, NANGLS) 
FORMAT {5E 10.3) 


CARD 4: V1, V2, DV 
FORMAT (3F 10.3) 


CARD 5: IRPT 
FORMAT {I5) 


Definitions of these quantities will be discussed in Section 3.2. 


3.2 Basic Instructions 


The various quantities to be specified on each of the five control cards along with the fourteen 
optional cards (summarized in Section 3.1) will be discussed in this section. 


3.2.1 CARD 1: MODEL, ITYPE, IEMSCT, IMULT, M1, M2, M3, M4. Md, 
M6, MDEF, IM, NOPRT, TBOUND, SALB 
FORMAT (1315, F8.3, F7.2) 


MODEL selects one of the six geographical-seasonal model atmospheres or specifies that user- 
defined meteorological data are to be used. 
MODEL =O _ If meteorological data are specified (horizontal path only) 

1 Tropical Atmosphere 
2 Midlatitude Summer 
3 Midlatitude Winter 
4 Subarctic Summer ae 
6 Subarctic Winter 
6 1976 US Standard 
7 


If a new model atmosphere (e.g. radiosonde daia) is to be read in, 
(NOTE: MODEL = 0 Used for horizontal path only) 
ITYPE Indicates the type of atmospheric path. 


ITYPE =i Fora horizontal (constant-pressure) path 
2 Vertical or slant path between two allitudes 
3 Fora vertical or slant path to space 


IEMSCT Determines the mode of execution of the program. 


IEMSCT=0 _ Program execution in transmittance mode 
1 Program execution in thermal radiance mode 


2 Program execution in radiance mode with solar/lunar single scattered 
radiance included 


3 Program calculates directly transmitted solar irradiance 
IMULT Determines execution with multiple scattering 


IMULT=0O _ Prograin executed without multiple scattering 
1 Program executed with multiple scattering 


(NOTE: IEMSCT must equal 1 or 2 for muitiple scattering) 


M1, M2, M3, M4, MS, and M6 are used to modify or supplement the altitude profiles of temperature and 


pressure, water vapor, ozone, methane, nitrous oxide and carbon monoxide from the atmospheric 
models stored in the program. 


MDEF Uses the default (U.S. Standard) profiles for the remaining species. 


For normal operation of program (MODEL 1 to 6) 
Sel M1=M2=M3=0, M4=M5=M6=MDEFs0 


Thesc parameters are reset to default values by MODEL (1 to 6) when they are equal to zero 


When M1=0 M1 reset to ‘MODEL’ 
When M2=0 M2 reset to 'MODEL' 
When M3=0 M@& reset to ‘MODEL’ 
When M4=0 M4 reset to ‘MODEL’ 
When M5=0 Mb reset to ‘MODEL’ 
When M6=0 M6 reset to 'MODEL' 
When MDEF=0 MDEF reset to 1 for all remaining species (not needed with MODEL I to 6). 


lf MODEL=0 or 7 and if: 


a. M1 through M6 are zero then the JCHAR parameter on card 2C.1 should be utilized to 
supply the necessary amounts. 
or 
b. M1 through M6 are non-zero then the chosen default profiles will be utilized provided the 
specific JCHAR option is blank. 


Ml=1to6 Default temperature and pressure to specified model atmosphere 

M2=1 to6 Default HgO to specified model atmosphere 

M3=1to6 Default Og to specilied model atmosphere 

M4=1to6 Default CH, to specified model atmosphere 

M5=1to6 Default NgO to specified model atmosphere 

M6=1to6 Default CO to specified model atmosphere 

MDEF=1 Use default profiles for COg, Og, NO, SO2g, NOg, NH3, HNOg, (not needed 
with MODEL 1 to 6). 


If MODEL=0 or MODEL=7, the program expects to read user supplied atmospheric profiles. Set: IM=1 
for first run. To then rerun the same user-aimosphere for a series of cases set IM=0; LOWTRAN will 
reuse the previously read data. 


IM=O For normal operation of program or when subsequent calculations are to 
be run with the MODEL data set lasi rcad in. 
1 When user input data are to be read initially. 
NOPRT=0_ For normal operation of program. Contrels TAPE6 output 
1 To minimize printing of transmittance or radiance table and 
atmospheric profiles 
-1 Controls TAPE 8 output 
TBOUND =Boundary temperature (°K), used in the radiation mode (if IEMSCT=1 or 2) 
for slant paths that intersect the earth or terminate at a grey boundary 
(for example, cloud, target). If TBOUND is left blank and the path 
intersects the Earth, the program will use the temperature of the first 
atmospheric level as the boundary temperature. 
SALB  =Surface albedo of the Earth at the location and average frequency of the 
calculation (0.0 to 1.0). If SALB is left blank the program assumes the 
surface is a blackbody (with enissivily equal to 1; for example, SALB=Q). 
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Table 7 summarizes the use of the five control parameters MODEL, ITYPE, IEMSCT, IMULT, and 
NOPRT on CARD 1. 


Table 7. LOWTRAN CARD 1 Input Parameters: MODEL, ITYPE, IEMSCT, IMULT and NOPRT 


A Pt 


MODEL, ITYPE, IEMSCT, IMULT, M1, M2, M3, 
CARD |] M4, M5, M6, MDEF, IM, NOPRT, TBOUND, SALB 
FORMAT (1315, F8.3, F7.2) 


COL COL COL COL COL *: 
a MODEL ITYPE IEMSCT IMULT NOPRT . 
Without : 
User- Horizontal Transmittance multiple 1 | TAPE 8 
defined* scattering Output 
. With TAPE 6 
Tropical Slant path Radiance : multiple Normal 
H1 to H2 scattering output 


Radiance with TAPE 6 | 
9 Midlatitude Slant path solar/lunar 1 Short 
summer to space scattering output : 
3 Midlatitude } 3 Transmitted 
winter 5 solar irradiance 


Subarctic 
summer 


Subarctic 
winter 


1976 U.S. 
Standard 


User- 
defined* 


M1, M2, M3, M4, M5, M6, MDEF, IM, TBOUND, SALB are left blank for standard cases. 
*Options for non-standard models. 


; 3.2.2 CARD 2: IHAZE, ISEASN, IVULCN, ICSTL, ICLD, IVSA, VIS, WSS, 
WHH, RAINRT, GNDALT 
FORMAT (615, 5F10.3) 


IHAZE, ISEASN, IVULCN, and VIS select the altitude and seasonal-dependent aerosol profiles 
and aerosol extinction coefficients. IHAZE specifies the aerosol model used for the boundary-layer 
(0 to 2 km) and a default-surface meteorclogical range. The relative humidity dependence of the 
boundary-layer aerosol extinction coefficients is based on the water vapor content of the model 
atmosphere selected by MODEL. ISEASN selects the seasonal dependence of the profiles for both the 
tropospheric (2 to 10 km) and stratospheric (10 to 30 km) aerosols. IVULCN is used to select both the 
profile and extinction type for thie stratospheric aerosols and to determine transition profiles above 
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the stratosphere to 100 km. VIS, the meteorological range, when specified, will supersede the default 
meteorological range in the boundary-layer aeresol profile scl by IHAZE, 

IHAZE selects the type of extinclion and a default meteorological range for the boundary-layer 
aerosol models only. If VIS is also specified, it will override the default IHAZE value. Interpolation of 
the extinction coefficients based on relative humidity is performed only for the RURAL, MARITIME, 
URBAN, and TROPOSPHERIC coefficients used im the boundary layer (O to 2 kin altitude). 

IHAZE=0 | no aerosol attenuation included in the calculation 

=l RURAL extinction, default VIS = 23 lan 

=2 RURAL extinction, default VIS = 5 km 

=3 NAVY MARITIME extinction, sets own VIS (wind and relative humidity dependent) 

=4 MARITIME extinction, default VIS = 23 km (LOWTRAN model) 

=5 URBAN e ‘inction, default VIS = 5 km 

=6 TROPOSPHERIC extinction, default VIS = 50 km 

=7 User defined aerosol extinction coefficients. Triggers reading cards 2D, 2D1 and 2D2 
for up to 4 altitude regions cf user defined exlinclion, absorption and asymmetry 
parameters 


=8 FOGI1 (Advective Fog) extinction, 0.2-km VIS 
=9 FOG2 (Radiative Fog) extinction, 0.5-km VIS 
=10 DESERT extinction sets own visibility from wind speed (WSS) 
ISEASN selects the appropriate seasonal aerosol profile for both the tropospheric and 


stratospheric aerosols. Only the tropospheric aerosol extinction coelflicients are used with the 2 to 10 
km profiles, 


ISEASN=0 _ season determined by the value of MODEL; 
SPRING-SUMMER for MODEL = 0, 1, 2, 4, 6, 7 
FALL-WINTER for MODEL = 3, 5 
=l1 SPRING-SUMMER 
=2 FALL-WINTER 
The parameter IVULCN controis both ihe selection of the acrosol profile as well as the type of 
extinction for the stratospheric aerosols. It also selects appropriate transition profiles above the 
stratosphere to 100 km. Meteoric dust extincifom coefficients are always used for altitudes from 30 to 
100 km. 
IVULCN=0,1 BACKGROUND STRATOSPHERIC profile and extinction 
=2, MODERATE VOLCANIC profile and AGED VOLCANIC extinction 
=3 HIGH VOLCANIC profile and FRESH VOLCANIC extinction 
=4 HIGH VOLCANIC profile and AGED VOLCANIC extinclion 
=5 MODERATE VOLCANIC profile and FRESH VOLCANIC extinction 
=6 MODERATE VOLCANIC profile and BACKGROUND STRATOSPHERIC extinction 
=7 HIGH VOLCANIC profile and BACKGROUND STRATOSPHERIC extinction 
=8 EXTREME VOLCANIC profile and FRESH VOLCANIC extinction 


Table 8 shows the value of IVULCN corresponding to the different choices of extinction coefficient 
model and the vertical distribution profile. 


Table 8. LOWTRAN CARD 2 Input Parameter: IVULCN 


VERTICAL DISTRIBUTION 


BACKGROUND MODERATE HIGH EXTREME 
STRATOSPHERIC VOLCANIC VOLCANIC VOLCANIC 


BACKGROUND 
STRATOSPHERIC 


AGED 
VOLCANIC 


EXTINCTION MODEL 


FRESH 
VOLCANIC 


ICSTL is the air mass character (1 to 10}, only used with the Navy maritime model (IHAZE = 3). 
Default value is 3. 


ICSTL=1 open ocean 


10 strong continental influence 


4 


ICLD specifies the cloud models and rain models used. 
The rain profiles decrease linearly from the ground to the top of the associated cloud model. 
The program cuts off the rain at the cloud top. 
ICLD=0 No clouds or rain 
=1 Cumulus cloud; base 0.66 km, top 3.0 km 
=2 Altostratus cloud; base 2.4 km, top 3.0 km 
=3 Stratus cloud; base 0.33 km, top 1.0 km 
=4 Stratus/Strato Cu; base 0.66 km, top 2.0 km 
=5  Wimbostratus cloud; base 0.16 km, top 0.66 km 


=6 2.0 mm/hr Drizzle (modeled with cloud 3) 
rain 2.0 mm/hr at O km to 0.22 mm/hr at 1.5 km 
=7 5.0 mm/hr Light rain (modeled with cloud 5) 
rain 5.0 nun/hr at O km to 0.2 mm/hr at 2,0 km 
=S 12.5 mm/hr Moderate rain (modeled with cloud 5) 
rain 12.5 mm/hr at 0 km to 0.2 mm/hr at 2.0 km 
=9 25.0 mm/hr Heavy rain (modeled with cloud 1) 
rain 25.0 mm/hr at 0 km to 0.2 mm/hr at 3.0 km 
=10 75.0mm/hr Extreme rain (modeled with cloud 1) 
rain 75.0 mm/hr at 0 km ta 0.2 mm/hr at 3.5 km 
=11 Read in user defined cloud extinction and absorption. Triggers reading Cards 2D, 
2D1 and 2D2 for up to 4 altitude regions of user defined extinction, absorption, and 
asymmetry parameters 
=18 Standard Cirrus model 
=19 Sub-visual Cirrus model 
=20 NOAA Cirrus model (LOWTRAN 6 Model) 


IVSA selects the use of the Army Vertical Structure Algorithm (VSA) for aerosols in the boundary 
layer. 


IVSA=0 not used 
=] Vertical structure algozithm 


VIS specifies the surface meteorological range” 36:37.38 (km) overriding the default value associated 
with the boundary layer chosen by IHAZE. If set to zero uses default value specified by IHAZE. 


* The terms "meteorological range” and "visibility" are not always used correctly in the literature. 
Correctly,36.37 visibility is the greatest distance at which it is just possible to see and identify 
with the unaided eye: (a) in the daytime, a dark object against the horizon sky; and (b) at night, 
a known moderaiely intense light Source. Meteorological range is defined quantitiatively, 
eliminating the subjective nature of the observer and the distinction between day and night. 
Meteorological range V is defined by the Koschmieder formula 


1, 1 3912 
Veen === 
Bove” B 


where B is the total extinction cocflicient which is the sum of the molecular and aerosol 
exlinction, and e is the threshold contrast, set equal to 0.02. As used in the LOWTRAN 
computer code, the inputs are in terms of meteorological range with B, the extinction 
cocflicient, evaluated at 0.55 pm, If only an cbserver visibility Vobs is available, the 
meteorological range can be estimated as V = (1.3 £0.3) - Vopg. (See Gordon?) 


36. Huschke, R-E. (1959) Glossary of Metecrology, American Meteorological Society, Boston, MA, 
p. 638 


37. Middleton, W.E.K. (1952) Vision Through the Atmosphere, University of Toronio Press, p. 250. 
38. Gordon, J.I. (1970) Daytime Visibility, A Conceptual Review, AFGL-TR-79-0257, AD AOS5451. 


VIS > 0, user specified surface meteorological range (km) 
= 0, uses the default meteorological range set by IHAZE (See Table 10). 


WSS specifies the current wind speed for use with the Navy maritime and desert aerosol models. 


WSS= current wind speed (m/s). Used with the Navy maritime model (IHAZI:=3) or the 
DESERT model (IHAZE=10). 


WHH specifies the 24 hour average wind speed for use with the Navy maritime model. 
WHH= 24-h average wind speed (m/s). Only used with the Navy maritime model (IHAZE=3) 


For the Navy Maritime mode! if WSS=WHH=O, default wind speeds are set according to the value of 
MODEL, see Table 9. For the Desert aerosol model (IHAZE=10), if WSS<0, the default wind speed is 
10 m/s. 


Table 9. Default Wind Speeds for Different Model Atmospheres Used with the 
Navy Maritime Model (IHAZE=3) 


WSS and WHH 
Default Wind Speed (m/s) 


User-defined (Horizontal Path) 
Tropical 

Midlatitude summer 
Midlatitude winter 

Subarctic summer 

Subarctic winter 

U.S. Standard 


User-defined 


RAINRT Specifies the rain rate 


RAINRT = Rain rate (mm/hy) default value is zero. 
Used to top of cloud when cloud is present, 
when no clouds rain rate used to 6km 


GNDALT specifies the altitude of the surface relative to sea level 


GNDALT = Altitude of surface relative to sea level (km) 
Used to modify aerosol profiles below 6 km altitudes 


Table 10 summarizes the use of the control parameters IHAZE, ISEASN and IVULCN on card 2 and 
Table 11 summarizes the use of the parameter ICLD. 


27 


Table 10, LOWTRAN CARD 2 Input Parameters: IHAZE, ISEASN, [VULCN, VIS 


IHAZE, ISEASN, IVULCN, ICSTL, ICLD, 
IVSA, VIS, WSS, WHH, RAINRT, GNDALT 
FORMAT (615, 5F 10.3) 


ISEASN IVULCN 
COL COL PROFILE/ 
ION} 10 |SEASON] 15 |SEASON| PROFILE JEXTINCTION |EXTINCTION 


NO AEROSOLS 
Set by Set by : 
model model 
Spring- Spring- Meteoric |’ 

2 5 1 Jsummer summer dust 
extinction 
Fall- Fall- 
3 2 winter winter 


Navy 
maritime 


[eo Background | Background Normal 
strato- strato- atmospheric 
maritime spheric spheric profile 
Tropospheric 
profile / Moderate Aged 
tropospheric volcanic volcanic 


extinction 


Fresh 
voleanic 


High 
voleanic 


Transition 


profiles 
High Aged - volcanic 
voleanic volcanic to nomnal 


Mocerate Fresh 
volcanic voleanic 
Moderaie 
voleanic 


Background 
strato- 
spheric 


Background 
High strato- 
volcanic spheric 


Fresh 
volcanic 


Extreme 
voleanic 


——e 


30 to 100 km 


Oto2 km 10 to 30 kin 


2 to 10 km 


* Default VIS, can be overridden by VIS > 0 on CARD 2 
** Sets own default VIS 
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fable 11. LOWTRAN CARD 2 Input Parameter: ICLD 


a 


IHAZE, ISEASN, IVULCN, ICSTL, ICLD, IVSA, 
VIS, WSS, WHH, RAINRT, GNDALT 
FORMAT (615, 5F 10.3) 


F 


OR CLOUD AND OR RAIN OPTION 


NO CLOUDS OR RAIN 
CUMULUS CLOUD 
ALTOSTRATUS CLOUD 
STRATUS CLOUD 
STRATUS/STRATO CUMULUS 
NIMBOSTRATUS CLOUD 
2,0 MM/HR DRIZZLE (MODELED WITH CLOUD 3) 
5.0 MM/HR LIGHT RAIN (MODELED WITH CLOUD 5) 
12.5 MM/HR MODERATE RAIN (MODELED WITH CLOUD 5) 
25,0 MM/HR HEAVY RAIN (MODELED WITH CLOUD 1) 
75.0 MM/HR EXTREME RAIN (MODELED WITH CLOUD 1) 
USER DEFINED CLOUD EXTINCTION AND ABSORPTION 
STANDARD CIRRUS MODEL 
SUB VISUAL CIRRUS MODEL 
NOAA CIRRUS MODEL (LOWTRAN 6 MODEL) 


3.2.2.1 Optional Cards Following CARD 2 


Optional input cards after CARD 2 selected by the parameters ICLD, IVSA, MODEL, and IHAZE on 
CARDS 1 and 2, 


CARD 2A CTHIK, CALT, CEXT, ISEED (if ICLD=18, 19 or 20) 
FORMAT {3F 10.3, 110) 

Input card for cirrus altitude profile subroutine when ICLD = 18, 19, or 20 

CTHIK is the cirrus thickness (km) 


If CTHIK=0 use thickness statistics 
>O user defined thickness 


CALT is the cirrus base altitude (km) 
CALT=0_ use calculated value 
>O user defined base altitude 


CEXT is the extinction coefficient (km7!) at 0.55 pm 


CEXT=0_ use 0.14* CTHIK 
>O user defined extinction cocfficient 


ISEED Is the random number initializing flag. 
ISEED=0, use default mean values for cirrus. 
>O, initial value of seed for random number generator, function RANF (SEED), 
(different values of SEED produce different random number sequences). This 
provides for statistical determination of cirrus base altitude (CALT) and thickness 
(CTHIK). 


NOTE: Random number generator is system dependent, 


CARD 2B ZCVSA, ZTVSA, ZINVSA (IfIVSA = 1) 
FORMAT (3F 10.3) 


Input card for Army VSA subroutine when IVSA = 1. The case is determined by the parameters 
VIS, ZCVSA, ZTVSA, and ZINVSA. 


CASE 1: cloud/fog at the surface; increasing extinction with height from cloud/fog base to 
cloud/fog top. Selected by VIS ¢ 0.5 km and ZCVSA 2 0. 


Use case 2 or 2' below cloud and case 1 inside it. 


CASE 2: hazy/light fog; increasing extinction with height up to the cloudbase. Selected by 
0.5 < VIS <10 km, ZCVSA 2 0. 


CASE 2': clear/hazy; increasing extinction with height, but less so than case 2, up to the 
cloudbase. Selected by VIS » 10 km, ZCVSA 2 0, 

CASE 3: no cloud ceiling bui a radiation fog or au inversion of boundary layer present; 
decreasing extinction with height up to the height of the fog or layer. Selected by 
ZCVSA < 0, ZINVSA 2 0. 


CASE 4: no cloud ceiling or inversion layer; constant extinction with height. Selected by 
ZCVSA < 0 and ZINVSA < 0. 
ZCVSA is the cloud ceiling height (km): 
i ZCVSA > 0.0 the known cloud ceiling height; 
= 0,0 height unknown: the program will calculate one for case 2, and 
default is 1.8 kim for case 2'; or 
< 0.0 no cloud ceiling (cases 3 and 4). 
ZIVSA is the thickness of the cloud (case 2) or the thickness of the fog at the surface ‘ 
(case 1) (km): 
If ZTVSA > 0.0 the known value of the cloud thickness; 
= 0.0 thickness unknown: default is 0,2 km, 
ZINVSA is the height of the inversion or boundary layer (km): 
If ZINVSA > 0.0 the Known height of the inversion layer; 
= 0.0 height unknown: default is 2 km, 0.2 km for log: 
< 0,0 no inversion layer (case 4, if ZCVSA < 0.0 also). 
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OPTIONAL USER INPUT CARDS 2C, 2C1, 2C2 and 2C3 


The following cards handle user input data. 
Cards 2C and 2C1 are always read for MODEL 0 or 7. 


CARD 2C ML, IRD1, IRD2, TITLE (MODEL=0/7, IM=1) 
FORMAT (315, 18A4) 


Additional atmospheric model (MODEL 0/7) 


New model atmospheric data can be inserted provided the parameters ‘MODEL’ and ‘IM’ 
are set equal to 0/7 and 1 respectively on card 1. 

ML = Number of atmospheric levels to be inserted (Maximum of 34) 

IRD1 Controls reading WN20, WCO... and WNH3, WHNOS3 (CARD 2C2) 


IRD1=0 No read 
‘ IRD1=1 Read CARD 2C2 


IRDN2 Controls reading AHAZE, EQLWCZ, ... (CARD 2C3) 
IRD2=0 No read 
IRD2=1 Read CARD 2C3 
TITLE = Identification of new modei atmosphere 
CARD 2C1 ZMDL,P,T, WMOL({1), WMOL(2}), WMOL(3), (JCHAR(J), J=1, 14) 
FORMAT (F 10.3, 510.3, 1841) 


CARD 2C2. (WMOL\), J=4, 12) 


FORMAT (8E 10.3) 
ZMDL = Altitude of layer boundary (im) 
P = Pressure of layer boundary 
T = Temperature of layer boundary 
WMOL (1-12) = Individual molecular species (see Table 11A. for species) 
JCHAR(1-14) = Control variable on units selection for profile inpui (P,T and molecular 


constituents, see Table 11A.) 

By utilizing a choice of values for the JCHAR(J) control variable (where J = 1,14) the user can 
designate specific uniis or accept defaulls for the various molecular species and for the temperature 
and pressure. If JCHAR(J) is left blank the program will default to the values chosen by M1, M2, M3, 
M4, M5, MG and MDEF when the given amount is zero. 


For JCHAR(1) 


A indicates Pressure in (mb) e 
B indicates Pressure in (atm) 
Cc indicates Pressure in (torr) 


1-6 will default to specified atmospheric value 
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Table 114. The Association of the JCHAR (J) Subscript (J=1,14) 
with the Variables P,T and WMOL 


VARIABLE SPECIES 


; 


P pressure 

T temperature 
WMOL{1) water vapor (IIgO) 
WMOL{2) carbon dioxide (COg) 
WMOLK{3) ozone (G3) 
WMOL{4) nitrous oxide (N9O) 
WMOL{5) carbon monoxide (CO) 
WMOL{6) methane (CH4q) 
WMOL{7} oxygen (Oo) 
WMOL{8) nitric oxide (NO) 
WMOL{9) sulphur dioxide (SO) 
WMO L{10) nitrogen dioxide (NO) 
WMOL(1 1) ammonia (NH3) 


WMOL(12) nitric acid (HNO3) 
— 


OON A AF WH 


ee ee ol ool 
PON K © 


For JCHAR(2) 


A indicates Ambient temperature in deg (K) 
B indicates Ambient temperature in deg (C) 
1-6 will default to specified atmospheric value 


For JCHAR(3- 14) 


A indicates Volume mixing ratio (ppmv) 

B indicates Number density (cm~°) 

C indicates Mass mixing ratio (gm/kg) 

n indicaies Mass density (gm/m°) 

E indicates Partial pressure (rab) 

F indicates Dew point temp (TD i T (K)) - HO only 

G indicates Dew point temp (TD in T (C)) - HgO only 

H indicates Relative humidity (RH in percent) - H20 only 
I is available for user definition 

1-6 will default to specified model atmosphere 


CARD 2C3 is read when IRD2 is set to 1 on CARD 2C. 


CARD 2C3 AHAZE, EQLWCZ, RRATZ, IHA1, ICLD1, IVUL1, ISEAIL, ICHR1 
FORMAT (10X, 3F 10.3, 515) 


AHAZE Aerosol or cloud scaling factor (equal to the visible [wavelength of 0.55 ym] 
extinction coefficient [km7!] at altitude ZMDL) 
[NOTE: only one of AHAZE or EQLWCZ is allowed] 
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EQLWCZ Equivalent liquid water content (gm/m?%) at allitude ZMDL for the aerosol, cloud or 
fog models 
RRATZ Rain rate (mm/h) at allitude ZMDL 


Only one of IHA1, ICLD1 or [VULI is allowed 


THAI Aerosol model extinction and meteorological range control for the allitude, ZMDL, 
See IHAZE (CARD 2) for options 


ICLD1 Cloud extinction control for the altitude, ZMDL, See ICLD (Card 2) for options. 
When using ICLD1 it is necessary to set ICLD (CARD 2) to the same value as the initial 
input of ICLD1. 


IVUL1 Stratospheric aerosol profile and extinction control for the allitude ZMDL, see 
IVULCN (CARD 2) for options 


The precedent order of these parameters (IHA1, ICLD1 and IVUL]) is as follows: 


If {IHA1>0) then others ignored 
If (HA1=0) and (ICLD1>0) then use ICLD1 
If (("HA1=0) and (ICLD1=0) then use IVUL1 


Wf AHAZE and EQLWCZ arc both zero, default profile loaded from IHA1, ICLD1, IVUL1 
ISEA1 Aerosol season control for the altitude, ZMDL, see ISEASN (CARD 2) for options, 


ICHRI Used to indicate a boundary change beiween 2 or more adjacent user defined aerosol or 
cleud regions at. altitude ZMDL (required for IHAZE=7 or ICLD=1}). 
ICHR1=0, no boundary change in user defined aerosol or cloud regions (regions are not 
adjacent). 
=1, signifies the boundary change in adjacent user defined aerosol or cloud 
regions. 


NOTE: ICHR1 internally defaulis to O if (IHA1¥7) or (ICLD1¥#1]1). 
OPITONAL CARDS 2D, 2D1, 2D2 


The foliowing cards allow the user to specify their own attenuation coefficients for any or all 
four of the aerosol regions. They are only read if IHAZE=7 or ICLD=11 are specified on card 2 (pages 24 
and 25}. 


CARD 2D (IREG (Ij, W=1,4) (QHAZE=7 or ICLD=11 input) 
FORMAT (415) 


IREG Specifies which of the four altitude regions a user defined aerosol or cloud model is used 
(IHAZE=7 /ICLD=11) 


(NOTE: Regions default to 0-2, 3-10, 11-80, 35-100 km and can be overridden with IHA1' settings in 
MODEL 7) 


IREG (Il) =O Use defauil values for region II 
IREG (1) = 1 Read extinction, absorption, and asymmetry for a region 


CARD 2D1 AWCCON, TITLE 
FORMAT (E10.3,18A4) 

AWCCON is a conversion factor from equivalent liquid water content (zin/m4) to extinction 
coefficient (km~!), It is numericaliy equal to the equivalent liquid water content 
corresponding to an extinction coefficient of 1.0 km7!, at a wavelength of 0.55 pm. 
AWCCON has units of (kin - gm - ni), 

TITLE for an aerosol or cloud region (up to 72 characters) 


CARD 2D2 (VX(I), EXTC(N, I), ABSC(N, I), ASYMIN, 1), 1=1,47) 
FORMAT (3(F6.2, 2F7.5, F6.4) ) 


Where N:-il when IREG(ID=1 for up to 4 altitude regioris. User defined aerosal or cloud extinction and 
absorption coefficients when IHAZE=7 or ICLD=11 


VX(I) = Wavelengths for the aerosol or cloud coefficients (not used by program: bul should be the 
same as the wavelengths defined in array VX2 in subroutine EXTDTA, see Table 12) 


EXTC(N, 1) = Aerosol or cloud extinction coefficierits, normalized so that EXTC for a 
wavelength of 0.55 pm (I=4) is 1.0 km“), 
ABSC(N, 1) = Aerosol or clu.id absorption coefficient, normalized so that EXTC for a 


wavelength of 0.55 pm (I=4) is 1.0 km7}. 
ASYM(N, I) = Aerosol or cloud asymmetry parameter 
3.2.3 CARD 3: H1, H2, ANGLE, RANGE, BETA, RO, LEN 
FORMAT (6F10.3, 15) 


CARD 3 ts uscd to define the geometrical path parameters for a given problem. 


Hi] = initial altitude (km) 
H2 = final altitude (km) (for ITYPE = 2) 
H2 = tangent height (km) (for ITYPE = 3) 


It is important to emphasize here that in the radiance mode of program execution (IEMSCT = 1 or 2), 
H1, the initial altitude, always defines the position of the observer (or sensor). H1 and 112 cannot be 
used interchangeably as in the transmittance mode. 


ANGLE = initial zenith angle (degrees) as measured from H1 

RANGE = path length (km) 

BETA = earth center angle subtended by H1 and H2 (degrees) 

RO = radius of the earth (km) at the particular latitude at which the calculation is to 
be performed. 
W RO is left blank, the program will use the midlatitude value of 6371.23 km if 
MODEL is set equal to 7. Otherwise, the earth radius for the appropriate 
standard model atmosphere (specified by MODEL) will be used as shown in 
Table 13. 


For an ITYPE = 2 path for which 111 > H2 (and by necessity, ANGLE > 90°), two paths are possible: the 
long path from I¥1 through a tangent height to H2 and the short path from H1 to H2. LEN selects the 
type of path in these cases. 


Tabie 12. The VX Array with the Required Wavelengths for the 


Multiply Read Card 2D2 
1 .2000 25 9.0000 
2 .3000 26 9.2000 
3 10.0000 
4 10.5910 
5 11.0000 
6 11.5000 
7 12.5000 
8 14.8000 
9 15.0000 
10 16.4000 
11 17.2000 
12 18.5000 
13 21.3000 
14 25.6000 
15 30.0000 
16 40.0000 
17 50.0000 
18 42 60.0000 
19 43 80.0000 
20 44 100.0000 
21 45 150.0000 
22 46 200.0000 
23 300.0000 


i] 
Ast 


Table 13. Default Values of the Earth Radius for Different Model Atmospheres 


Earth Radius, RO (km) 


Model Atmosphere 


Not used 
6,378.39 


User-defined (Horizontal Path) 


Tropical 


| Midlatitude summer 6,371.23 
Midlatitude winter 6,371.23 
Subarclic summer 6,356.91 
Subarctic winter 6,356.91 
U. S. Standard 6,371.23 


User-defined 6,371.23 


LEN = 9 © short path (default). 


1 long path through the tangent height. 


il 


It is not necessary to specify every variable on CARD 3; only those that adequately describe the 
problem according to the parameter ITYPE (as described below). (See Table 14). 


(1) Horizontal Paths (ITYPE = 1) 
(a) Specify H1, RANGE 


(b) if non-standard meteorological data are to be used, that is, if MODEL = O on CARD 1, then 
refer to Section 3.3 for a detailed explanation. 


(2) Slant Paths Between Two Altitudes (ITYPE = 2) 
(a) specify H1, H2, and ANGLE 
{b) specify H1, ANGLE, and RANGE 
(c) specify H1, H2, and RANGE 
(d) specify H1, H2, and BETA 
(3) Slant Paths to Space (ITYPE = 3) 
{a) specify H1 and ANGLE 


(b) specify H1 and H2 (for limb-viewing problem where H2 is the tangent height cr 
minimum altitude of the path trajectory). 


For case 2(b), the program will calculate H2, essuming no refraction; then proceed as for case 2(a). 
The actual slant path range will differ from the input value. This method of defining the problem 
should be used when refraction effects are not important; for example, for ranges of a few tens of km at 
zenith angles less than 80°. For case 2(c), the program will calculate BETA and then proceed as for case 
2(d). For case 2(d), the program will determine the proper value of ANGLE (including the effects of 
refraction) through an iterative procedure, This method can be used when the geometrical 
configuration of the source and receiver is known accurately, but the initlal zenith angle is not known 
precisely due to atmospheric refraction effects. Beta is most frequently determined by the user from 
ground range information. 


Table 14 lists the options on CARD 3 provided to the user for the different types of atmospheric 
paths. 


3.2.3.1 Alternate CARD 3 for Transmitted Solar or Lunar Irradiance ; 
(IEMSCT = 3) 


For calculating directly transmitted solar or lunar irradiance, an ITYPE = 3 path is assumed and 
CARD 3 has the following form: 

ALTERNATE CARD 3 

H1, H2, ANGLE, IDAY, RO, ISOURC, ANGLEM 

FORMAT (3F 10.3, 15, 5X, F10.3, 15, F10.3) 


altitude of the observer a 
tangent height of path to sun or moon 
apparent solar or lunar zenith angle at H1 


Hl 


i 


as 
z 
a 
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Table 14. Allowable Combinations of Slant Path Parameters 


a ee Se 


RANGE BETA 


LEN option is available only when Hi > H2 and ANGLE > 90°. Otherwise, LEN is 
set in the program, 


(2) H2 calculated assuming no refraction. Calculated RANGE will differ from the 
input value. 


(3) BETA calculated assuming no refraction. 
(4) Exact ANGLE is calculated by iteration of the path calculation. 


(5) H2 is interpreted as the tangent height. if HZ and ANGLE are both zero, Case 3a is 
assumed with ANGLE = 0 (that is vertical path). For a path tangent at the earth's 
surface, read in a small number for H2, for example, 0.001 km. 


IDAY = day of the year, used to correct for variation in the earth-to-sun 
distance 

RO = radius of earth (default according to MODEL) 

ISOURC = O extraterrestrial source is the sun 


= 1 extraterrestrial source is the moon 
ANGLEM = phase angle of the moon, that is, the angle formed by ithe sun, moon, 
and earth (required if ISOURC = 1) 


Either H2 or ANGLE should be specified. If both are given as zero, then a vertical path (ANGLE = 
0°) is assumed. If IDAY is not specified, then the mean earth to sun distance is assumed. 

If the apparent solar zenith angle is not known for a particular case, then the solar scattering 
option (IEMSCT = 2) may be used along with, for instance, the observers location, day of the year, and 
time of day to determine the solar zenith angle (see section 3.2.3.2 of the user instructions). Note that 
the agparent solar zenith angie is zenith angle at H1 of the refracted path to the sun and is less than 
the astronomical solar zenith angle. The difference between the two angles is negligible for angles less 
than 80°, 


3.2.3.2 Optional Cards Following CARD 3 


Optional input cards after CARD 3 are selected by parameters IEMSCT on CARD 1 and IPH on 
CARD 3A1. 
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CARD3A1_TPARM, IPH, IDAY, ISOURC (if IEMSCT = 2) 
FORMAT (415) 


Input card for solar/lunar scattered radiation when IEMSCT = 2. 


IPARM 
IPH 


IDAY 


ISOURC 


0,1,2 controls the method of specifying the solar/lunar geometry on CARD 3A2. 

O Henyey-Greenstein aerosol phase function (see CARD 3A2) 

1 uSer-supplied aerosol phase function (see CARD 3B) 

2 MIE-generated internal database of aerosol phase functions for the LOWTRAN 

models 

day of the year, that is, from 1 to 365 used to specify the earth to sun distance and 
(if IPARM = 1) to specify the sun's location in the sky. (Default value is the mean 

earth to sun distance, IDAY=93). 

O extraterrestrial source is the sun 

1 extraterrestrial source is the moon 


CARD 3A2 PARM1, PARM2, PARM3, PARM4, TIME, PSIPO, ANGLEM, G 
(IEMSCT = 2) 
FORMAT (8F 10.3) 


Input card for solar/lunar scattered radiation when IEMSCT = 2, Definitions of PARM1], 
PARM2, PARM3, PARM4 determined by value of IPARM on CARD 3A1, (See Table 14A,). 


PARM 1 


PARM2 
PARM3 
PARM4 


NOTE: 


PARM 1 
PARM2 
PARM3, 


For IPARM = 0 


observer latitude (-90° to +90°) 

(Note that if ABS(PARM 1) is greater than 89.5° the observer is assumed to be at 
either the north or south pole. In this case the path azimuih is undefined. The 
direction of line-of-sight must be specified as the longitude along which the path 
lies. This quantity rather than the usual azimuth is read in.) 


observer longitude (0° to 360°, west of Greenwich) 
source (sun or moon) latitude 


source (sun or moon) longitude 


For IPARM = 1 


The parameters IDAY and TIME must be specified. This option cannot be used with 
ISOURC = 1. 


observer latitude ( -90° to +90°) 
observer longitude (0° to 360°, west of Greenwich) 


PARM4 are not required 
(Note that the calculated apparent solar zenith angle is the zenith angle at H1 of the 
refracted path to the sun and is less than the astronomical solar zenith angle. The 
difference between the two angles is negligible for angles less than 80 degrees.) 


Table 14A. Card 3A2; Options for Different Choices of IPARM 


Observer Latitude 
(-90 to +90°) 


Observer Latitude 
(-90 to +90°) 


Azimuth Angle 
Between Observer LOS 

& Observer to Sun Path 
Solar Zenith Angle 


Observer Longitude Observer Longitude 


(O to 360° West of (O to 360° West of 
PARM2 Greenwich) Greenwich) 


Source Latitude 


Source Longitude 


Greenwich Time 
(Decimal Hours) 
Path Azimuth Angle 
(degrees East of Due 

North) 


Path Azimuth Angle 
(degrees East of Due 
North) 


ANGLEM 
(only if ISOURC=1) 


Lunar Phase Angle Luaar Phase Angle 


| Asymmetry Asymmetry Asymmetry 
G Parameter Parameter Parameter 
(only if IPH=0O) (-1 to +1) (-1 to +1) (-1 to +1) 
for use with Henyey- for use with Henyey- for use with Henyey- 
Greenstein Phase Greenstein Phase Greenstein Phase 
Function Function Function 
For IPARM = 2 
PARM1 = azimuthal angle between the observers line-of-sight and the observer-to-sun path, 


measured from the line of sight, positive east of north, between -180° and 180° 
PARM2 = the sun's zenith angle 
PARM3, PARM4 are not required 


REMAINING CONTROL PARAMETERS 


TIME = Greenwich time in decimal hours, that is, 8:45 am is 8.75, 5:20 pm is 17.33 etc. 
(used with IPARM = 1) 


PSIPO = path azimuth (degrees east of north, that is, due north is 0,0”, due east is 90.0” etc. 

(used with IPARM = 0 or 1) 

phase angle of the moon, that is, the angle formed by the sun, moon, and earth 

(required only if ISOURC = 1) 

G = asymmetry factor for use with Henyey-Greensteiu phase function (only used with 
IPH = 0), e.g., +1 for complete forward scattering, O for isotropic or symmetric 
scattering, and -1 for complete backscaitering. 


B 
© 
Cc 
ie} 
K4 
" 


CARD 3B1 NANGLS (Only if IPH = 1 on card 3A1) 
FORMAT (15) 


Input card for user-defined phase functions when IPH = 1. 
NANGLS = _ number of angies for the user-defined phase functions (maximum of 50) 


CARD 3B2 (1 to NANGLS) 
(ANGF(), F(1,0, F(2,0), F(3.1), F(4,.0), I= 1, NANGLS) 
FORMAT (5E 10,3) 


Input card for user-defined phase functions when IPH = 1. 
ANGFI(I) = scattering angle in decimal degrees (0.0° to 180.0°) 


F(1,0 = user-defined phase function at ANGF(I), boundary layer (0 to 2 km default 
altitude region) 

F(2,)) = user-defined phase function at ANGF{I), troposphere (2 to 10 km default 
altitude region) 

F(3,)) = user-defined phase function at ANGF()), siralusphere (10 io 30 kin default 
altitude region) 

F(4,)) = User-defined phase function at ANGF(I), mesosphere (30 to 100 km default 


altitude region) 
The default altitude regions may be overriden by the parameters IHA1, ICLD1 or IVUL1. 


3.2.4 CARD 4: V1, V2, DV 
FORMAT (3F 10.3) 


The spectral range and increment of the calculation. 
Vi = initial frequency in wavenumber (cm!) 
V2 = final frequency in wavenumber (cm7!), where V2 > V1 


DV = frequency increment (or step size) (em!) 
(Note that v = 104/4, where v is the frequency in cm! and ) is the wavele agth in 
uum, and that DV can only take values that are a multiple of 5 cn! and that V1 and 
V2 are reset to the next lowest or highest integer mulliple of 5 cm!) 


3.2.5 CARD 5: IRPT 
FORMAT (15) 


The control parameter IRPT causes the program to recycle, so that a series of problems can be 
run with one submission of LOWTRAN., 


IRPT = Oto end program 
= 1 toread all new data cards (1, 2, 3, 4, 5) 
= 2not used 
= 3 read new CARD 3 (the geometry card) and CARD 5 
= 4read new CARD 4 (frequency) and CARD 5 
4 or IRPT = 2 will cause program to STOP 


Thus, if for the same model atmosphere and type of atmospheric path the reader wishes to 
make further transmittance calculations in different spectral intervals V1' to V2’ etc., and fora 
different step size (DV etc.), then IRPT is set equal to 4. In this case, the card sequence is as follows and 
can be repeated as many times as required. 


CARD5 IRPT=4 
CARD 6 VI'V2'DV' 
CARD 7 IRPT=4 
CARD & V1" V2" DV" 
CARD 9 IRPT=0 


The final IRPT card should always be a blank or zero. When using the IRPT option, the 
wavelength dependence of the refractive index is not changed (use the IRPT = 1 option if this is 
required), 

Table 15 summarizes the user-control parameters on CARD 4 and CARD 5. 


Table 15. LOWTRAN CARD 4 and CARD 5 Input Parameters: V1, V2, DV, iRPT 


CARD 4 V1, V2, DV 
Format (3F 10.3) 
P| V1 (cm!) V2 (em!) DV (cm™}) Multiple of 5 cm™! 


CARD 5 IRPT 
Format (15) 


IRPT 
End of program. 
Read new CARDS 1, 2, 3, 4, and 5. 
Not used (stops program). 
Read new CARDS 3 and 5. 
Read new CARDS 4 and 5. 


| cos | 
0 
1 
2 
3 
4 


NOTE:  IRPT=3 cannot be used when running multiple scattering cases or 
solar single scattering. Use IRPT=1. 


3.3 Non-Standard Conditions 


Several options and combination of choices are available if atmospheric transmiltance/ 
radiance caiculations are required for non-standard conditions. Here non-standard refers to 
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conditions other than those specified by the parameters MODEL, IHAZE, and ICLD on CARDS 1 and 2. 
These options enable the user to insert: 


(1) An additional atmospheric model (MODEL = 7), which can be in the form of radiosonde or 
other source data. It is not necessary to duplicate the altitudes used in LOWTRAN 7, 
(2) Meteorological conditions for a given horizontal path calculation (MODEL = 0) ; o 
(3) A combination of any or all of the 12 gases can be input for cach layer boundary with defailt : 
choices interleaved with user supplied data. 
(4) Aerosol vertical distributions can be input at specified altitudes by Unc use of AHIAZE, 
EQLWC2Z, and/or IHAi on CARD 2CS when IRD2 ts set to 1 on CARD 2C. 
(5) Cloud liquid water contents and or rain rates can be input at specified altitudes by the use of 
EQLWCZ, RRATZ and/or ICLD 1 on CARD 2C3 when IRD2 is set to 1 on CARD 2C. 


(6) Any combination of the one to four Aerosol altitude regions can be replaced by reading in 
specific values of extinction and absorption coefficients and asymmetry parameters for specific 
regions by utilizing CARDS 2D, 2D1 and 2D2. The parameters can be for aerosols and for clouds. 


3.3.1 ADDITIONAL ATMOSPHERIC MODEL (MODEL = 7) 


A new model atmosphere can be inserted by the use of CARD 2C and the required multiples of 
card 2C), provided the paraineters MODEL and IM are set to 7 and 1 respectively on CARD 1, The 
number of atmospheric levels to be inserted (ML) must also be specified 011 CARD 2C. anee 


The appropriate meteorological parameters and the format are given below: 


CARD 2C ML, IRD1, IRD2, TITLE 
FORMAT (315, 18A4) 


CARD 2C1 ZMDL, P, T, WMOL(1), WMOL{2), WMOL{3), (JCHAR(W), J = 1,14) 
FORMAT (F 10.3, 5E10.3, 15A1) 


See section 3.2.2.1 above for a detailed description of each variable. 


3.3.2 HORIZONTAL PATHS (MODEL = Q) 


If known meteorological data are to be used for horizontal path atmospheric transmittance / 
radiance calculations, then set MODEL = 0 and IM = 1 on CARD 1. Proceed to read the meteorological 
conditions utilizing CARDS 2C and 2C1 as described above. In this instance the parameter ML must be : 
set to 1. 


3.3.3 USER INSERTED VALUES FOR ATMOSPHERIC GASES 
(MODEL = O OR 7) 


The user may wish to enter specific values of any or all of the atmospheric gases. This can be 
accomplished by utilizing CARDS 2C, 2C1 and 2C2. Ou CARD 2C set IRD1 = 1. The specific gas 
amounts for individual gases can be entered on CARDS 2C! and 2C2 and by utilizing the parameter 
JCHAR on CARD 2C1, the user has a choice of entering data in several dillerent sets of units (eg, 


volume mixing ratio, number density, .. . etc), or defaulting to one of the model atmospheric gases al 
the specified altitude. 


3.3.4 USER INSERLED VALUES FOR AEROSOL VERTICAL DISTRIBUTION 
(MODEL = 0 OR 7) 


The capability exists for the user to be able to replace aerosol distributions at specific allitudes. 
In order to accomplish this the user must set IRD2 to 1 on CARD 2C. Then specify the altitudes on 
CARDS 2C1 along with the variables, defaults and or units by utilizing the parameters as explained in 
section 3.2.2.1 

On CARD 2C3 the aerosol scaling factor for a given allilude can be entered by using the variable 
AHAZE, or an appropriate value for EQLWCZ, or defaulted by using the variable IHAL. 


3.3.5 USER INSERTED VALUES FOR CLOUD AND OR RAIN RATES 
(MODEL = 0 OR 7) 


The same capability exists permitting the user to replace cloud liquid water contents and or rain 
rates at specific altitudes as described in the above section. This is accomplished by setting IRD2 to 1 
on CARD 2C. Then the specific altitudes may be entered on CARDS 2C1 along with the variables, 
defaults and or unit selection by using the remaining parameters of CARD 2C] as described earlier. 

On CARD 2C38 the variables EQLWCZ and/or RRATZ can be used to enter the intended value of 
equivalent liquid walter content of a cloud and/or the rain rate at the specified altitude, or the cloud 
attenuation can be specified by using AHAZE. The uscr may default at the specilied altitude to one of 
the built-in cloud and/or rain model vai'es by using ICLD1. 


3.3.6 REPLACEMENT OF AEROSOL OR CLOUD ATTENUATION MODELS 
(IHAZE = 7 AND/OR ICLD = 1]) 


The aerosols or cloud model utilized in any or all of the four altitude regions may be replaced by a 
user input model. The built-in regions are 0-2 km, > 2-10 km, > 10-30 km and > 30-100 kin, These 
regions may be modified by the use of the parameters IHA1, ICLD1i or IVULI, This option is initialized 
by setting INAZE = 7 or ICLD = 11. 

On CARD 2D the variable IREG (1, 4) determines which of the altitude regions will have 
replacement values read in. The user is required to enter a conversion factor, AWCCON (km - gm - m-3), 
on card 2D1, which converts aerosol or cloud profiles specified in terms of equivalent Hquid water 
content, EQLWCZ (gm - m~9), to an extinction coefficient (km~!), This conversion factor (AWCCON), is 
only used if the aerosol or cloud concentration are specified by EQLWCZ instead of by the visible 
extinction, AHAZE, The LOWTRAN values for this variable are stored as DATA statements in 
subroutine EXABIN, ( See DATA statements ELWCR, ELWCM, ELWCU, ELWCT, AFLWC, RFLWC, 
CULWC, ASLWC, STLWC, SCLWC, SNLWC, BSLWC, FVLWC, AVLWC, and MDLWC.) 

The muluiply read CARDS 2D2 (13 cards) consist of four variables, VX, EXTC, ABSC and ASYM, 
The first variable VX is the wavelength of the data points which should correspond to the wavelengths 
used in the program (defined in array VX2 in Subroutine EXTDTA, see table 12), The next three 
variables EXTC, ABSC, and ASYM are the aerosol or cloud extinction, absorption coefficients and the 
asymmetry parameters respectively. As stated previously the variable IREG (1-4) will determine if (he 


uSer is reading in 1, 2, 3, or 4 sets of CARDS 2D1-2D2. Additionally, by utilizing the variables IVUL1 
and ISEAI the user can substilute for stratospheric acrosol profiles and can change the seasonal 
profile valucs. 

The values of EXTC(N,]) and ABSC(N,I) should be normalized so that EXTC(N,4) = 1.6 (Le., the 
extinction for wavelength 0.55 um is normalized to 1.0). 


4. EXAMPLES OF PROGRAM OUTPUT 


Four test cases are explained in this section. The input cards for the four cases are listed in Table 
16 and the listings of the file QUTPUT (TAPE 6) are in Tables 17 through 20. The next four subsections 
contain detailed explanations for the specific test cases, The output of case 1 is discussed in detail 
while the discussion of the next three cases is limited to a description of the physically important ‘ 
details of the input. The LOWTRAN 7 tape as distributed by the National Climatic Data Center (see 
Section 2.4) contains additional test cases. 


4.1 Case 1: Single and Multiple Scattered Solar Radiance 


The program will compute atmospheric radiance using the single scatlered solar radiance 
calculation combined with the multiple scattered solar radiance calculation. The parameters selecied 
for this case are as follows: the atmospheric profile is the 1976 US Standard and the boundary layer 
aerosol model is RURAL with 5 km meteorological range. Since MJ] through M6 are lett blank along 
with MDEF being set to O, the atmospheric profile valucs (pressure and temperature) along with all of 
the molecular species amounts are retained as initially chosen by MODEL, The surface albedo is set io 
0.4, The path for this case is from 20 km to 0 km looking straight down. The solar zenith angle is 60 
degrees and the relative azimuth angle is 0. The earth to sun distance is chosen for day 1 (January 1). 
The aerosol phase functions are from the Mie-generated database stored in the program. The spectral 
range of the calculation ts from 14000 to 34000 cn! in steps of 100 em™!. (This test case is admittedly 
under sampled and is for illustration only. If more spectral detail exists in a given region, a smaller 
DV is recommended). 

The output for this case shown in Table 17 will now be described in detail, On the first page the 
output echoes the inpui cards exactly as they are read in. CARD 1 selects the US Standard Atmosphere 
(MODEL = 6), type of slant path (ITYPE = 2), single scattering option (IEMSCT = 2), 1.ultiple scattering 
option (IMULT = 1), and the surface albedo (SALB = 0.4). CARD 2 selects the RURAL aerosol profile 
with the default meteorological range of 5 km (IHAZE = 2). CARD 3 describes the slant path in terms of 
the observer altitude (H1 = 20.0), endpoint at the ground (H2 = 0.0) and zenith angle (ANGLE = 180). On 
CARD 3A] IPARM = 2 determines that the sun position will be described on CARD 3A2 in terns of the 
relative azimuth at H1 of 0° (PARM1) and the solar zenith angle at H1 of 60° (PARM2). IPH = 2 specifies 
that the Mie-generated phase function will be used while IDAY = 1 specifies that the earth-to-sun 
distance used to calculate the incident solar intensity will be that for January 1. ISOURC = 0 specifies 
the sun as the source. On CARD 3A2, only PAIRM1 and PARM2, as described earlier, are used 
for this case. The remaining parameters are left blank, CARD 4 selecis the spectral range of 14000 to 
34000 cm! in steps of 100 cem-!. Following CAD 4 the program interprets the input cards in a more 
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readable fashion. 

The next three pages labeled Atmospheric Profiles, list the profiles of pressure, temperature, and 
absorber densities at the indicated altitudes Z for the selected atmospheric model. The units of the 
densities are as indicated or are in relalive units (-), The densities of the nitrogen continuum. (Ng) are 
in scaled LOWTRAN units. The self-broadened water vapor continuum density factor ‘CNTMSLF' is 
listed next, The column labeled 'MOL SCAT is the total air density relative to 1013.25 mb and 273,15K 
which is used to calculate molecular scattering. 'N-1' is the refractivity or the index of refraction -1.0. 
The ultra-violet ozone and oxygen [Og (UV), Og (UV)] densities are the aciual densilies. 

The second of these three pages contains the following information: the foreign-broadened 
water vapor continuum density factor 'CNTMFRN' and the nitric acid 'HNO3' density factor. AEROSOL 
1 to 4 refer to the aerosol scaling factors for the 4 aerosol altitude regimes. "AER1*RH' is the log 
weighted relative numidity weighted by the boundary layer aerosol. 'CIRRUS' is the cirrus cloud 
density profile in km7!, Since the cirrus cloud option was not selected, the cirrus density is zero 
everywhere. The last column labeled 'RH' is the path averaged relative humidity. 

The third page of these profiles contains the densities for each molecular species in the following 
order? HgO, Og, COg, CO, CHy, NgO, Og, NH3, NO, NOo, and SOg. The next six tables labelled: Slant Path 
Parameters in Standard Form, Calculation of the Refracted Path Through the Atmosphere, 
Cumulative Absorber Amounts for the Path from H1 to Z, Summary of the Geometry Calculation, 
Equivalent Sea Level Total Absorber Amounts and Single Scattering Point to Source Paths are 
repeated twice in the output. The first set of tables result from the calcuiations taken for a vertical 
path through the entire atmosphere (O to 100 km). This calculation is needed to obtain the up and 
down fluxes necessary to perform multiple scattcring calculations. The sccond set of tables are for the 
line of sight path, in this case 20 to O km, 

The first table contains pertinent information concerning the slant path parameters. All 
possible combinations of slant path parameters are reduced to the standard set of H1, H2, ANGLE, PHI, 
HMIN and LEN, where PHI is the arrival angle at H2, HMIN its the altitude of the lowest point along the 
path and LEN has a value of unity if the path goes through a tangent height and a value of zero 
otherwise. 

A layer-by-layer description of the ray trace through the almosphere appears next. The ray tract 
always begins at the lowest point along the path HMIN and proceeds upwards. The various quantities 
presented are: layer boundaries; zenith angle at the bottom of the layer THETA; curved path length 
through the layer DRANGE; cumulative path length RANGE :rom HMIN; earth-centered angle 
subtended by the layer DBETA and by the cumulalive path BETA; arrival angle at the top of the layer 
PHI, path bending within the layer DBEND and cumulative BENDING: densily-weighted path averaged 
pressure and temperature PBAR and TBAR for the layer; and average total air density RHOBAR. 

The next table lists the cumulative absorber amounts for the species described in the 
atmospheric profile table from the observer to the different layer boundaries. 'CNTMSLF1' and 
‘CNTMSLF2' refer to the self-broadened water vapor amount and (he temperature correction factor. A 
sunnary of the geometry calculation appears next, followed by “Equivalent Sca Level Total Absorber 
Amounis" that are the total amounts for the whole path for each of the aforcementioned variables, 
"Mean RII" refers to the mean relative humidity in the boundary layer only. 


Next is the output of the single solar scattering geometry routines: this includes a sununary of 
the path from cach scattering point to the sun. ‘SUBTENDED ANGLE’ is the earth-ccntered angle 


between the observer at H1 and the scattering point. The ‘SOLAR ZENITH’ angle is the astronomical or 
unrefracted zenith angie to the sun, However, if the refractive bending for the path to the sun is 
calculated to be greater than 0.1°, then the solar zenith angle is corrected by the bending and that path 
is calculated again (up to four times) until the correction for refraction is less than 0.1°. 

The 'RELATIVE AZIMUTH’ is negative because internally the azimuth angles are measured 
positive counterclockwise, while for the input cards, azimuth angles are specified posilive clockwise. 
'SCTTR ANGLE' is the scattering angle at that point. A message will be printed if any scattering point 
is in the shade. The last column contains the molecular phase function used for the calculation. 

AS previously noted the same six tables are now repeated, only this time the calculations are 
performed for the chosen path, in this case 20.0 to 0.0 km. 

Finally the results of the calculation are printed, In this case, it is the radiance and the total 
transmittance for the path, The first two columns are the frequency and wavelength in wavenumbers 
(cm!) and micrometers (4m) respectively. The "Total Radiance" is broken down into three types: 
"Atmospheric Radiance”, which includes all thermal radiation emitted by the atmosphere or by the 
boundary (and thermal radiation scattered by the atmosphere or reflected by the ground); "Path 
Scattcred Radiance”, which includes the solar radiation scattered by the atmosphere (the single 
scattered component of the solar radiation, "S Scat", is also shown separately); "Ground Reflected 
Radiance", which includes the solar radiation reflected by the ground (both the directly transmitted 
solar reaching the ground, "Direct", which is shown separately, or the solar radiation scattered within 
the atmosphere before reaching the ground). The “Total Radiance" is the sum of the three types of 
radiances reaching the observer. The four major components: atmospheric radiance, path scattered 
radiance, ground reflected radiance and total radiance are listed per cm~! and um. The next column 
‘INTEGRAL’ is the integrated radiance from the initial frequency. The last column lists the total 
transmittance from H1 to H2. A short summary is printed at the end of the table, including the actual 
boundary temperature and emissivity used (if any), the integrated absorption, average transmittance, 
integrated radiance and the minimum and maximum values of the total radiance, which are useful in 
determining the limits for plotting. 

The last ilem printed is the value of IRPT on CARD 5. In this case IRPT = O indicates there are no 
further input cards following. Figure 6 is a depiction of this case, the solid line is the total multiple 
scattered radiance and the dotted line shows the contribution due to single scattered solar radiance. 


4.2 Case 2: Thermal Radiance 


This case which is shown in Figure 7 consists of two separate runs of LOWTRAN, The solid line is 
the result using the muitiple scattering calculations and the dotted line is the result of a standard 
radiance calculation. The output is depicted in Table 18. 

In both runs MODEL has been set to 2 choosing the midlatitude summer model atmosphere, 
ITYPE = 3 is a Slant path to space, IEMSCT = 1 is the thermal radiance option and M1 through M6 are 
left blank and MDEF set to O retains the amounts for each gas as chosen by MODEL. The only change 
on CARD 1 is the IMULT parameter which is on for one run and off for the second. On CARD 2 IHAZE 
has been set to 2 thereby choosing the rural aerosol profile with a meteorological range of 5 km. In 
defining the geometrical path Hl and H2 have been sel to 0, while angle has been set to 80°. Since this 
is a slant path to space the program will perform calculations between O to 100 km, HMIN will be set to 


(Text continued on page 93) 
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O and the RANGE will be calculated (for this case 479.887 km). The spectral range for these 
calculations have been chosen as 740 to 1250 cm”! with steps of 5 en}, 


4,3 Case 3: Atmospheric Transmittance 


‘This case represents a standard atnicspheric transmittance run with the added feature of having 
an elevated surface of 1.5 km. The parameters selected for this case are as follows: the atmospheric 
profile is the 1976 US Standard and the boundary layer aerosol model is RURAL with 23 km 
meteorological range. Since M1 through M6 are left blank along with MDFEF being set to O, the 
atmospheric profile (pressure and temperature), along with all of the remaining molecular species are 
retained as initially chosen by MODEL. By setting the variable GNDALT to 1.5 an elevated surface is 
chosen. The path is a vertical path from the ground to 6 km. The spectral range of the calculation is 
from 990 to 1090 cm! in steps of 5 cm7). Figure 8 is a plot from 750 - 1250 cn}, containing tne results 
from this case. The TAPE 6 output appears in Table 19. 


4.4 Case 4: Directly Transmitted Solar Irradiance 


Case 4 calculates the directly transmitted solar irradiance from 6000 to 7500 cm*! for an 
observer at the ground and a solar zenith angle of 60°. This case is selected by IEMSCT = 3 on CARD 1. 
An ITYPE = 3 slant path is required: CARD 3 specifies the observer altitude H1 = 0.0, the apparent 
solar zenith angle at H1, ANGLE = 60°, and the day of the year IDAY = 74 (the Ides of March). The TAPE 
6 output appears on Table 20. 

The spectral output consists of the transmitted solar irradiance and the incident solar 
irradiance, corrected for variations in the earth-to-sun distance according to the day of the year. Also 
given are the cumulative integrals of these two quantities and the total transnuttance for the path 
from H1 tothe sun. The directly transmitted solar irradiance (solid line} and the solar radiance 
(dashed line) is depicted graphically in Figure 9. 


4.5 Description of TAPE7 and TAPE8 


Two files are opened for the storage of pertinent output, UNIT = 7 and UNIT -. 8. The next section 
will describe the various types of output included on TAPE 7. 


4.5.1 TAPE7 OUTPUT 


The spectral results from each successful case run are written to a file on UNIT = 7, along with a 
11-line header identifying the parameters for the case. This file can then be used to further proce.ss the 
results. For example, the plot. filter, and scanning function programs described in Appendices A, B, 
and C work from this file. 


The data on the file consists of an 11-line header followed by the spectral data. The header 
consists for the most part of images of the input control cards, with the default values supplied for 
parameters that have defaults, for example VIS on CARD 2. The first four lines of the header are 
CARD 1, CARD 2, CARD 2A, and CARD 2B, When the cirrus, VSA, and single scattering options are net 


(Text continued on page 114) 
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selected, the images of their respective control cards are supplied with the parametcrs set equal to -99, 
The fifth line of the header is the image of CARD 2C and gives the number of atmospheric layers and 
title of the atmospheric profile. This information is given for all values of MODEL (0 to 7) even though 
it is read in only for MODEL = 0 or 7, CARDS 2C1, 2C2, 2C3, 2D1 and 2D2 are not included in the 
header. Line 6 contains the values of the path parameters on CARD 3 after they have been evaluated. 
CARDS 3B1 and 3B2 are not included. Lines 7 and 8 are the images of CARDS 3A1 and 3A2, with all the 
parameters set equal to -99 if the single scattering option is not selected. Line 9 is the imag: of CARD 
4, the frequency card, and line 10 gives the value of IRPT from CARD 5. Line 11 consists of headers for 
the columns of data to follow. 

The following lines cont 1in the spectral cutput and differ depending on the type of calculation: 
transmiltance, radiance, radiance with single scattering, or directly transmitted solar irradiance. 
For a transmittance calculation (IEMSCT = 0) each line of spectral data has the following format: 


frequency (cm7!) F7.0 
transmittance, total 

transmittance, water vapor (band type) 

transmittance, uniformly mixed gases 

transmittance, ozone 


transmiltance, trace gases 
transmittance, nitrogen continuum 
transmittance, water vapor continuum 11F8.4 


OPN oan O DN & 


trausmiilance, molecular scattering 
transmittance, aerosol and hydrometeor 
transmittance, nitric acid 

aerosol and hydrometeor absorptance 
Log of the total transmittance 1PE10.3 


a 
ones 


For a thermal radiance run (IEMSCT = 1) the format is: 


1. frequency (em7}) F7.0 
2. transmittance, total F8.4 
3. atmospheric radiance (W cm? ster {em7!)-}) 1PE9.2 
4. Log of the total transmittance 196, £10.3 i 
For a radiance case with single scattering or multiple scattering (IEMSCT = 2) the format is: 

1. frequency (cm7!) F7.0 : 
2. transmittance, total ee ¥F8.4 
3. atmospheric radiance (W em~2 ster” }(em=!)-}) 
4. path scattered radiance (W cm-? ster~}(em=!)-}) 
5. single scattered radiance (W cm~ ster }{cm=!)7!) 1P6E9.2 bi 
6. total ground reflected radiance (W cm~ ster” (enr“})-3j _ 
7. direct reflected radiance (W cm~2 ster7(em7})7}) : 
8. total radiance (W cm”? ster™(em™!)~}) 
9, TEBI - OP2F8.4 

10. TEB2SV | 

ll. Log of the total transmittance See T96, 1PE10.3 


TEB] is the total trarismittance for the L-shaped path from H1 to H2 to the sun. TEB2SV 1s the total 


transrnittance for the path from H1 to the last boundary before H2 to the sun. 


For directly transmitted solar irradiance (IEMSCT = 3) the format is: 
1. frequency (cm7!) 

transmittance, total 

transmitted solar irradiance (W cm~? {em~!)-}) 

incident solar irradiance (W cm-2 (cm7})-}) 

Log of the total transmittance 


The end of the spectral data is marked by a frequency of -9999. 


ahkwp 


4.5.2 TAPE8 OUTPUT 


F7.0 

F8.4 
1PE9.2 
E9.2 

T96, F10.3 


There are three different types of spectral data sent to UNIT=8, depending on cases run. The first 
type, printing of the separate transmittances of the molecular absorbers, only occurs with IEMSCT = 
O. In this instance the spectral results will be printed as follows: the same 10 line header printed ta 
UNIT = 7 will be sent to UNIT = 8 as well. Two lines of descriptive headers follow. Then the spectral 


data will have the following format: 


— 


frequency (cm~!) 

transmittance, water vapor (band type) 
transmittance, ozone 
transmittance, carbon dioxide 
transmittance, carbon monoxide 
transmittance, methane 
transmittance, nitrous oxide 
transmittance, oxygen 
transinittance, ammonla 
transmittance, nitric oxide 
transmittance, nitrogen dioxide 
12. transmittance, sulphur dioxide 


OC ON pP &h 


st pe 
- © 


The end of the spectral data is marked by a frequency of -9999. 


F7.0 


11F8.4 


The second type, printing of the black body function and differential transmittance for each 


Ja, along the line of sight and each frequency, occurs when NOPRT = -1, IEMSCT = 


1 or 2, and 


IMULT =0. The spectral results are stored as follows: the same 10 line header sent to UNIT = 7 will be 
sent to UNIT = 8 also. Line 11 is a line of descriptive headers for the data following. The output 
consists of a line of data for each boundary and repeats this procedure for each spectral frequency. 


The spectral results have the following format: 


1. frequency (cem~!) 

2, first altitude of boundary (km) 

3. second altitude of boundary (km) 

4. black body function (W cm ster”! (um)7!) 
5. differential transmittance 

6. transmittance, total 


F10,0 
F7.2 
F7.2 
IPE12.5 
E£12,5 
F12.9 


eS 


The third type, printing of the upward and downward, solar and thermal vertical fluxes, and 
solar irradiances for each layer and each frequency occurs when NOPRT = -1, IEMSCT = 1 or 2 and 
IMULT = 1, Again the same ten line header sent to UNIT = 7 will be sent to UNIT = 8. The next line 
consists of descriptive headers for the data io follow. The output consists of a line of data for each 
layer and repeats this procedure for each spectral frequency. The spectral results have the following 


format: 
1. frequency (cm7!) F10.0 
2. altitude of layer (km) F7.2 
3. upward total flux (W cm? (um)-}) 
4. upward solar flux (W cm”? (um)-!) 
5. downward total flux (W cm? (um)~!) 1P5E12.5 ‘ 
6. downward solar flux (W em-2 (um)-!) 
7. direct solar irradiance (W cm~? (uum)7}) - 
8. transmittance, total F105 


The end of the spectral data is marked by a frequency of -9999. 

When using this method to output data to UNIT = 8, it is necessary to follow a specific procedure. 
The position of the sun should be placed directly behind the viewer's back. This is accomplished by 
choosing the IPARM = 2 option on CARD 3A1 and then setting PARM1 = 180.0 on CARD 3A2, The 
second requirement is that the sum of the two variables, ANGLE of line of sight (CARD 3) and the solar 
zenith angle (PARM2 on CARD 3A2), must be equaij to 186.0, These procedures are necessary in order 
to obtain the correct direct solar irradiance. 


11. 
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Appendix A 


LOWTRAN 7 Plot Program 


The LOWIRAN 7 ploi prograin, LOWPLT, is an independent program designed to plot the 
transmittance and/or radiance output of a LOWTRAN 7 run. The program will plot one file of data, 
several files of data, or the same file of data as many times as desired. Another feature of the plot 
program is the ability to choose any of the many variables produced by a transmittance run, radiance 
run, or a radiance with scattering run. Finally, the user is given the option of plotting these variables 
on individual plots or on the same plot. When plotting different variables on the same plot, it is 
important to define the proper scale for the Y-axis. The LOWTRAN 7 plot program uses Tape 7 or 
Tape 8 as generated by a LOWTRAN 7 transmittance run. 

This plot program utilizes plot routines unique to the CDC-Cyber CALCOMP plotter. Therefore, 
some minor changes may be necessary to accommodate a different computer plotting systern. (See 
Section A3 for more details on CALCOMP system differences). 


Al. INSTRUCTIONS FOR USING THE LOWTRAN 7 PLOT PROGRAM 


The plot program extracis the data to be plotted from Tape 7 or Tape 8. To initiate the progr 
five input cards are required. 


Al.1 Input Data and Formats 


The data necessary to specify a given plot are given by the following five cards: 


CARD 1 PROGID, SCALE 
(FORMAT (3A10, F10.4) ) 


CARD 2 XSIZE, PFRBEG, PFREND, DELTAX, ITYP, IXAXIS, NUMFIL 
(FORMAT (4F 10.4, 315) ) 
CARD 3 YSIZE, YRMIN, YRMAX, DELTAY, ICRV, IYAXIS, NMYDEC 
(FORMAT (F 10.4, 3E10.2, 315) ) 
CARD 4 IPOS 
(FORMAT (515) } 
CARD 5 ISAMFL, ISAMPT 


(FORMAT (215) ) 


Al1.2 Utilizing Input Cards 


To produce multiple plots; input cards 2 through 5 must be re, eated. To end plotting: the first 
value on the following input card 2 should be negative. 
Definitions of the input card variables will be discussed in Section A2. 


AZ. BASIC INSTRUCTIONS 


The various quantities lo be specified on each of the five control cards are defined in this section, 
If plotting the variables from Tape 8, see section A4. 


A2.1 CARD 1 PROGID, SCALE 


The variable PROGID is a 30 Hollerith character identification header that is printed as a banner 
at the start of the plot. 
SCALE = Multiplicative factor to increase or decrease the plot size (usually = 1.0) 


A2.2 CARD 2 XSIZE, PFRBEG, PFREND, DELTAX, ITYP, [XAXIS, NUMFIL 


XSIZE = Length of X-axis in inches 
PFRBEG = Beginning wavenumber on plot in em! or wavelength in jun 
PFREND = Ending wavenumber on plot in cm! or wavelength in jum 
DELTAX = The interval between the numbered values on the X-axis 
ITYP = 0 Radiance per pum vs pan a 
= ] Radiance per em! vs env! 
= 2 Transmittance vs wim 
an Transmittance vs em7! 
IXNAXIS = O X-axis will be linear 
= | X-axis will be logarithmic 
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NUMFIL = 


Uses next available file of data 
Uses file of data specified by NUMFIL 


42.3 CARD 3 YSIZE, YRMIN, YRMAX, DELTAY, ICRV, IYAXIS, NMYDEC 


YSIXE = 
YRMIN = 
YRMAX = 


Length of Y-axis in inches 
Minimum transmittance or radiance value to be plotted 
Maximum transmittance or radiance value to be plotted 


Ina log plot YRMIN and YRMAX are input as the values of the respective exponents. In a linear plot 
they are entered as the actual minimum and maximum value, 


DELTAY = 
ICRV = 
> 
> 
< 
WYAXIS = 
NMYDEC = 


Oak WHR OO 


A2.4 CARD 4 IPOS 


If IEMSCT 


IPOS 


Arh wWwNr OO 


ND 


The Interval between the numbered values on the Y-axis 


Normal plot of line 

Calls special plotting routine to plot dashed and dotted lines 
Solid line without symbols 

Dashed line without symbols 

Dotted line without symbols 

Alternating dashes and dots without symbols 

Alternating dashes and 2 dots without symbols 

Same as ICRV = 1 to 5 with symbols at every point 
Alternating dashes of different lengths 

Data points only using symbol number = | ICRV |, where symbol 
numiber refers to the computer-system plotting table 


Y-axis will be near 
Y-axis will be logarithmic 


Number of digits to the right of the decimal point on the Y-axis 


Plot of total transmittance 

Plot of HyO transmittance ‘ 
Plot of uniformly mixed gases transmittance 
Plot of O3 transmittance 

Plot of trace gases transmittance 

Plot of Ng continuum transmittance 

Plot of H,O0 continuum transmittance 

Plot of molecular scattering transmittance 
Plot of aerosol/hydromieteor transmittance 
Plot of HNQg transmittance 

Plot of aerosol/hydrometeor absorptance 
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If IEMSCT =1, 


IPOS =O 


if TEMSCT — =2, 


IPOS =O 


for IPOS = 1-6 


Plot of total transmittance 
Plot of atmospheric radiance 
(units chosen by variable ITYP on Card 2) 


Plot of total transmittance 

Plot of atmospheric radiance 
Plot of path scattered radiance 
Plot of single scattered radiance 
Plot of ground reflected radiance 
Plot of direct reflected radiance 
Plot of total radiance 


The units are chosen by ITYP. 


NOTE: When running without multiple scattering (IMULT = 0) the path scattered radiance is equal to 
the single scattered radiance and the ground refleeted radiance is equal to the direet radiance. 


If IEMSCT =3, 
IPOS =O 
=] 

=2 


Plot of total transmittance 
Plot of transmitted solar irradiance 
Plot of incident solar irradiance 


for IPOS 1 and 2 the unils are chosen by ITYP. 


A2.5 CARD 5 ISAMFL, ISAMPT 


ISAMFL and ISAMPT are set to handle the next set of data to be plotted. ISAMFL is used when 
the user wishes to plot a second, third, or ,.., variable from the same file of LOWPRAN 7 Tape 7 data. 


The user should be sure that (he range of the variable falls within the range as specified on the plot 


axis by CARD 2 and CARD 3. 


ISAMFL =0 
=] 
ISAMPT =0 
=] 


Normal advance to next file of data 
Rewind and stay on same file 


Normal advance to new plot 
Plot data on same physical plot 


A3. THE AFGL COMPUTER SYSTEM CALCOMP PLOT ROUTINES 


The program LOWPLT uses several AFGL computer-systein CALCOMP plot calls Ghat may diver 


from the system plotting routines available to the user. It is anticipated (hat with suitable adaptition 


of the CALCOMF plot calls including the calls to PLTID3, PLOT, ENDPLT, NUMBER, SYMBOL, and 
LINE plotting can be accomplished with minimal difficulty. 


The general functions of the AFGL system). CALCOMP plot routines are as follows: 


1. CALL PLTID 3 (PROGID, KMAX, YMAX, FACTOR) 


This subroutine must be the first routine called as it initializes the plot. 


PROGID = A30-character hollerith array used as an identifier on the plot header 
XMAX = Max dimension in X inches (X limit for entire plot) (real) 
YMAX = Max dimension in Y inches (Y limit for entire plot) (real) 


FACTOR 


Amultiplicative factor to change size of plotting (usually 1.0) in both the X and Y 
directions 


2. CALL ENDPLT - PHYSICAL END OF PROGRAM 


Subroutine ENDPLT must be the last plotUing subroutine called in all levels. ENDPLT willemply 
plotting buffers and write identifying information at the end of the page. A call to ENDPLT causes 
termination of the job. 


3. CALL LINE (X, Y, N, K, J, L, XMIN, DX, YMIN, DY, SYMSZE) 


Subroutine LINE produces a single line by connecting the points defined in the dimensioned 
variables X and Y. 


X = Array cf X values (real) 
Y = Array ot Y values (real) 

N = Number of points to be plotted (integer) 

K = Repeat cycle (usually K = 1) (integer) 


When K = 2 the first, third, fifth, etc., points will be plotted 

When K = 3 the first, fourth, seventh, etc., points will be plotted, ete. 

Example CALL LINE (X(1), X(2), N2....) This is usually used when X and Y arrays are 
one mixed array. 


J = Control for using symbols (integer) 
Jo = Owill produce a line plot without symbols 
J = 1 will produce a line plot with a symbol at every point x 
J = 2will produce a line plot with a symbol at every second point, ete., a negative 
J will suppress the lines between the points. , 
L = Anumbez describing the symbol to be used. Only symbols whose integer equivalent 


is Oto 13 are centered around the point (X,Y). For all others the point (X.Y) is at the 
lower left corner of the symbol. 


XMIN = Starting value of X-axis, in units of X (real) 
DX = The difference between the final value on X-axis and the starting value of X-axis 
divided by the length of X-axis in inches (real). 
YMIN = Starting value of Y-axis, in units of Y (real) 
DY = The difference between ihe final value on Y-axis and the starting value of 


Y-axis divided by the length of Y-axis in inches. 
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SYMSZE = Anumber defining the size of the symbol to be used at the point (X,Y) (real). (Doctoral 
size is 0.08 in.) 


4. CALL PLOT (X, Y, IC) 


Subroutine PLOT is used to move the pen and to redefine a new origin, 
X = X Coordinate, in inches (real) 
Y = Y Coordinate, in inches (real 
IC = IfIC =2, pen down as pen moves to (X, Y) (integer) 
WIC = 3, pen up as pen moves to X, Y. 
If IC =-2 or -3, a new origin is defined at X, Y. 
NOTE: the pen will move to location X,Y on page in all cases. XA and Y are defined with respect to 
' the previously defined origin. Page frame limits should be considered in all cases. PIC ¥ 4.2 or 3, an 
error message is printed. 


5. CALL NUMBER (X, Y, HGHT, FPN, THETA, N) 
Subroutine NUMBER will interpret and plot a real (floating point) number. 
X = X Coordinate of lower left-hand corner of first digit, in inches, relative to the current 
origin (real) 
Y = Y Coordinate of lower left-hand corner of first digit, in inches, relative to the current 
origin (real) 
HGHT = Height of numbers to be plotted, in inches (real) 
FPN = Number to be plotted (real) 
THETA = Orientation of the number with respect to the X-axis, counterclockwise in degrees (real) 
N = Number of digits after the decimal point (integer), N =—1 will suppress the decimal point 
6. CALL SYMBOL (X, Y, HGHT, BCD, THETA, N) 
Subroutine SYMBOL will draw a series of symbols as defined in (he symbol table in the 
CALCOMP instruction manual. 
X = Coordinate of the lower-left corner of the first character, in inches, relative to the 
current defined origin (real) 
3 Y = Coordinate of the lower-left corner of the first character, in inches, relative to the 
current origin (real) 
‘ HGHT = The height of the characters, in inches (real). For pen plots the width of each character 
will be equal to the height. 
BCD = This parameter and the last parameter in CALL to SYMBOL (called N) determine the 


type of annotation the routine produces. If BCD is the text to be used as annotation, 
usually BCD or A type format, the characters must be left-justified and contiguous ina 
single variable, in an array, or in a Hollerith literal. Parameter N must eontain the 
number of characters to be plotted. IM BCD ts a single character of text, the text must be 
right-justified and parameter N = O. 

THETA = The angular orientation with respect to the X axis counterclockwise, degrees (real). 


pes 
i) 
a 


N = ‘This parameter plus parameter BCD detennines type of lettering/symbols produced by 
routine SYMBOL, 
N > O - defines character count in array BCD, left-justilied. 
N = 0 - defines single characters to be plotied, right-justified, 


N < O (negative) - determines the condition of the pen in the move from its present position to the place 
where the symbol is to be produced. 

IN =-1, the pen is up during the move, after which a symbol is produced. lf N = -2 or less, the pen is 
down during the move, after which a suymbol is produced. 

Three of the CALCOMP plot calls from the AFGL plotting Mbrary used in LOWPLT, are not 
consistent with the standard CALCOMP versions. These calls and their possible modifications are as 
follows: 

1. CALL PLTIDS - should be replaced with CALL PLOTS (0, 0, LDEV) where LDIbV is the plot 
device designation: 

2. CALL ENDPLT - should be replaced with a plot completion call appropriate to the user's 
system; and 

3. CALL LINE - a program modification must be made so that the starting value end the sealing 
factor gumediately follow the data values in the X-array and the Y-array respectively. 


A4, PLOTTING OF TAPES VARIABLES 


When plotting any of the variables on a Tape 8 output created by a transmittance run, the card 
sequence as outlined above remain the same. The only change is in the definition of the control 
variable IPOS on CARD 4 as follows: 


IPOS 


O Plot of HyO transmittance 


= 1 Plot of Og transmittance 


iT 
i) 


Plot of COy transmittance 
= 3 Plot of CO transmittance 
= 4 Plot of CHg transimitiance 
Plot of NYO transmittance 
Plot of Og transmittance 
Plot of NH transmittance 
Plot of NO transmittance 
Plot of NOg transinittance 


Plot of SOyg transmittance 


NOTE: The plot programu does not handle Tape 8 type output not created by a transmittance run. 


A5. LOWTRAN PLOT: PROGRAM STRUCTURE 


The plot program for LOWTRAN 7 is a separate program package available with LOWTRAN 7, The 
plot code structure is shown in Figure Al. Descriptions of each subroutine are shown in Table Al. 


LOWPLT 


CSHDOT | DSHDT2 | 


Figure Al. Plot Pregram Structure 
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LOWPLT 


AXISL 
AXLOG 
PLTDTA 


DRAW 


DASH 2 
DOT 
DSHDOT 
DSHDT2 


Table Al. Descriptions of Plot Program Subroutines 


a 


Main driver program. Reads 4 control cards. 


Plots a labeled linear axts. 

Plots a labeled logarithmic axis. 

Reads LOWTRAN data from TAPE 7. Delernnines which 
variable to plot. 


Sets up data to be plotted. Determines which type of line 
to plot. 


Plots a dashed line. 

Plots a dotted line, 

Plots a line with alternating dashes and dots. 
Plots line with alternating dash and two dots. 


Appendix B 


LOWTRAN 7 Filter Function Program, LOWFIL 


Bl. INTRODUCTION 


The LOWTRAN 7 Filter Function Program, LOWFIL, is primarily designed to calculate the 
effective atmospheric transmittance that would be measured by a filtered transmissometer. It can 
also be used to calculate the effective radiance seen by a radiometer measuring the radiance emitted or 
scattered by the atmosphere. 

This program is wriilen as an independent program package for use with the LOWTRAN 7 code. It 
assumes that the basic mass storage file, TAPE7, written by LOWTRAN, is available. It was written to 
be fairly fexibie in its use, although it is recognized that users may wish to modify the code to tailor it 
more specifically for their own applications. 

The effective average atmospheric transmittance measured by a transmissometer looking at a 
blackbody source is given by: 


Sa Jew) F(v) Biv. Typ) dv (Bl) 
IFO) Blv, Typ) dv 
where 


tT (v) is the spectral transmittance of the atmosphere, 
F (v) is the total combined response function of the sensor, including the transmittance of any 


filters, the effect of the optical components and the detector response, 
B(v, Typ) is the spectral emissivity of the blackbody source, 
Tp is the temperaiure of the blackbody, and 


vis the frequency and the range of integration includes all the non-zero values of F{v). 
The atmospheric or scattered solar radiance measured by a rad: »meter is proportional to: 


R= fico) Fvydv, (B2} 


where I (v) is the spectral radiance. In general, when using Eq. (B2) to determine the radiance 
measured by a radiometer it would be necessary to include in F (v), the radiometer angular field-of - 
view, calibration constants, and all other frequency independent quantitics that can be neglected for 
Eq. (B1), since they cancel out. The effective normalized radiance is given by: 


R —Ji(vyF(v)dv 
lee = F Sop (B3) 


max max 


where F_,,,, = the maximum value of F(v). When used with radiance oulput, LOWFIL will calculate both 
Rand Ig. The units of R will depend on the units of F(v), the sensor response function; Ip will be in 
W/cm2-ster, 

The effective filter weighted average transmittances for the different atmospheric components 
used in LOWTRAN are also computed, (for example, the transmittance due to ozone, the water vapor 
continuum, aerosol, etc.}. It should be recognized, that while: 


™(v) = TH ol Y) TCOo aoe bee THINO3(Y) ’ (B4) 
in general: 
q ed THyot COQ" i s -THNOS ‘ (B95) 


where the average transmittance for each of the atmospheric components is found by evaluating an 
expression analogous to Eq. (B1). 

In addition to the eight average component transmittances, the water continuum and band-type 
absorption are combined into a total water transmiltance. 


Twater(Y) = Too Band(?) TH20 cont, (Y) (B) 
Also the other gaseous transmittances are combined: 
TGaslY) a Toco, V) »TozonclY) « Tno(V) » Tal, scat(V) . THNO.,(Y) (37) 


These combined transmittance components hive proved useful in developing simple analytic 
expressions for the effective transmittance measured by a system. 


T= Twater + tTGas + TAcrosol + (B8) 


where analytic expressions are derived separately for twater + Teas + ANA Trgrosol « 
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B2. COMPUTATIONAL DETAILS OF THE PROGRAM 


The program starts by reading in a set of filter functions, or systems response functions. There 
can be up to 15 different filters each specified at up to 80 different wavelengths (these limits can be 
modified by adjusting the dimension statements} in each set. 

The program is written so that the systems response functions or filter function F(v), can be 
input either as a function of wavelength (um) or as function of wavenumber (cm~!), If they are input as 
a function of wavelength, the wavelengths are converted to wavenumbers and the order of the points is 
reversed to be in ascending order in wavenumber. 

Next, a specified number of LOWTRAN cases are read in from TAPE7 one at a time. For each 
LOWTRAN case, the different filters are cycled over one-by-one. The subset of the LOWTRAN 

) wavenumbers, which bracket the values for which the filter is defined, are found for each filter in 
turn, Then the values of the filter function are interpolated to this subset of the LOWTRAN 
) wavenuinbers. 

Finally, the integrals in Eqs. (B1) through (B3) are evaluated using the trapezoidal rule. After a 
set of filters is used with the specified number of LOWTRAN cases, the next set of filter functions is read 
in, and repeated until the input is terminated. 


B3. USE OF THE FILTER FUNCTION PROGRAM 


B3.1 Ge ®ral Remarks 


The filter function program LOWFIL utilizes TAPE7, the standard mass storage output file as 
written by the LOWTRAN 7 computer code. The control cards for the filter are read from TAPE 5, 
which in the listing that follows is equated to the INPUT file. 

The program is written so that a set of up to 15 filter functions (system response functions) can be 
read in at once, and used with one or several consecutive LOWTRAN calculations. The filter functions 
can be input either as a function of wavelength (jm) or as a function of wavenumber (cm™!). The 
wavelengths (wavenumbers) should be in ascending order, but can be at unequal intervals. An option 
is included to rewind the LOWTRAN TAPE 7 output file, so more than 15 filter functions (system 
response functions) can be used or the same filters can be re-used with a different blackbody source 
temperature. It should be emphasized that the spectral range of the LOWTRAN calculations should 
include the full sensitive range of the systems response functions to be used. 


B3.2 Control Cards 


There are three (3) basic control cards: 


CARD 1] NF, NEW, IFT, TEMP, IPRINT, NLOW 
(FORMAT (215, F10 2, 215)) 


Repeat CARD 2 and CARD 3 "NF" times, where NF is specified on CARD 1. 


CARD 2 IDFIL, KODE, IFWV, NW . 
(FORMAT (2A10, 315) “4 
CARD 3 (WAVE (I), FF (1), 1 = 1, NW) 


CARD 3 is a free format with as many cards as needed for all the values of WAVE and FF. 


B3.2.1 CARD 1 NF, NEW, IFT, TEMP, IPRINT, NLOW 
NF indicates number of different filters 


NF > 0O_ Read in NF filters (NF < 15) 
= QO Use preceding filter for next LOWTRAN output 
<0 Stop filter program 


NEW is an option to reuse the previous LOWTRAN data set for the next set of filters. 


NEW = O No, Read next NLOW sets of LOWTRAN data t 
= 1 Yes, Rewind the LOWTRAN TAPE 7 output file 


IFT is an option to enter blackbody temperature of source 


IFT = O No blackbody 
= 1 Fold in blackbody emmissivity 


TEMP = blackhody source temperature in degrees Kelvin 
IPRINT controls the information printed 


IPRINT 210, Print LOWTRAN transmittances and results below 
5, Print filter function with blackbody function folded in 
O, Only print the filter weighted transmittances 


VW 


NLOW indicates number of LOWTRAN data sets to use with this set of filters 


B3.2.2 CARD 2 IDFIL, KODE, IFWV, NW 


If NF is set to zero on CARD 1, a new input set of CARDS 2 and 3 are not read. If NF is > zero then : 
"NF" input sets of CARD 2 and 3 are read. i: 


IDFIL = 20 Hollerith character identification for the given filter 
KODE = _ Filter identification number (5 digits) 


IFWV is an option to input the filter function either vs wavelength or vs wavenumber, 


IFWV = = O Wavelength 
= 1 Wavenumber 
NW = Number of wavelengths or wavenumbers for the filter (NW < 80) 


B3.2.3 CARD 3 (WAVE (I), FF (J), I = 1, NW) 


The variable NW on input card 2 gives the number of wavelength or (wavenumber) filter function 
pairs necessary to CARD 3 or multiples of CARD 3. These should be in ascending order by wavelengtn 
(or wavenumber). 
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WAVE 
FF 


Wavelength or wavenumber, depending of IFWV 


Corresponding filter function 


B4. STRUCTURE OF THE FILTER FUNCTION PROGRAM 


The filter function program is an independent pregram package available for use on LOWTRAN 7 
TAPE7 output. The filter function code structure is shown in Figure Bl. Descriptions of each 
subroutine are iisted in Table Bl. 


LOWFIL 


INTGRT 


WAVEN BRACK INTLOG 


BIGFIL INTRAD BLKBDY 


INTGRT 


Figure B1. Filter Program Structure 
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Table Bl. Description of LOWTRAN Filter Program Subroutines 


| LOWFIL - 
WAVEN - 
BRACK - 
INTLOG - 


BIGFIL - 


INTRAD - 


BLKBDY - 
INTGRT - 


COMBT - 


Main driver program. Reads control cards and TAPE7 output from 
LOWTRAN 


Changes a system response function vs wavelength (11m) to a response 
function vs wavenumber (cm7!) 


Finds the LOWTRAN wavenumbers, which bracket the system response 
function 


Interpolates a pair of vectors, F{I) = f(x;), and X(I) = Xj, i= 1,2,...N,toa 
new set of coordinates; FNEW (J) = {{x;) and XNEW (J) = xj.J =1,2,...M 


Routine to find the maximum value of the system response function 


Integrates the emitted or scattered radiance from LOWTRAN times the 
systems response function 


Weights the system response function by a blackbody radiance 


Finds the average value of the transmittance from LOWTRAN weighted 
by the systems response function 


Finds the weighted average transmittance due to water vapor (both band- 
type and continuum combined) and all the other gases combined 
together 
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Appendix C 


LOWTRAN 7 Scanning Function Program 


The LOWTRAN 7 scanning function program is an independent program supplied with the 
LOWTRAN 7 package. The prograrn is designed to convolve a triangular slit function with the 
transmittance, or radiance available on the TAPE7 output file created by a LOWTRAN 7 run. The fixed 
half width at half maximum can be input as either wavenumbers (cm~!) or wavelength (um). This 
program allows the user to generate spectra with a coarser resolution than the current LOWTRAN 
resolution of 20 cm7!. Use of the scanning function preserves the correct integrated intensities or 
transmittances and avoids the errors of undersampling the spectra. 


C2. INSTRUCTIONS FOR USING THE SCANNING FUNCTION 
Four tape files are initially opened. The scanning function program exiracis the data to be 


scanned from TAPE7. TAPES is needed to read the control card, TAPE6, is used for printed output 
from the program and TAPE9 is used to store the scanned variable and is available for plotting. 


C2.1 Input Data and Format 


The data necessary to pass a scanning function over the TAPE7 variable is given by the following 
card: 


CARD 1 V10, V20, HWHM, IHW, NUMFIL, IVAR 
FORMAT (3F 10.5, 3:5) 


C3. BASIC INSTRUCTIONS 


The various quantities to be specified on the control card are defined in this section. 


C3.1 CARD 1 V10, V20, HWHM, IHW, NUMFIL, IVAR 


V10 = First output frequency in cm”! 

v20 = Last output frequency in cm”! 

HWHM~ =_ Fixed half width at half max (cm~!) defined at V10 (if IHW = 0) \ 
Fixed half width at half max (jum) (if IHW = 1) 

IHW = Flag for scanning function half width (0 or 1 as explained above) ‘ 

NUMFIL = Ouses next file of data 


> Ouses file of data specified by NUMFIL 


IVAR: The use of the variable VAR depends on the value of the variable IEMSCT as read from 
CARD 1 of TAPE7. 


If IEMSCT = 0 


IVAR = scans total transmittance 

scans H2O transmittance 

scans uniformly mixed gases transmittance 
scans Og transmittance 

scans trace gases transmittance 

scans No continuum transmittance 

scans HO continuum transmittance 

scans molecular scattering transmittance 


scans aerosol/hydrometeor transmittance 


i 
OMAN AA SP WN Ke © 


scans HNO; transmittance 
=10 scans aerosol/hydrometeor absorptance ‘ 


If IEMSCT = 1 


IVAR = 0 scan of total trarismittance 
1 


scan of atmospheric radiance 


es 


If IEMSCT = 2 


IVAR = O © scan of total transmiltance 
= 1 scan of atmospheric radiance 
= 2 scan of path scattered radiance 
= 3 scan of single scattered radiance 
= 4 © scan of total ground reflected radiance 
= 5 — scan of direct reflected radiance 
= 6 — scan of total radiance 


NOTE: When running without multiple scattering (IMULT = 0), the path scattered radiance is equal to 
the single scattered radiance and the ground reflected radiance is equal to the direct radiance. 


If lEMSCT = 3 
IVAR = scan of total transmittance 


.@) 
= 1 scan of transmitted solar irradiance 
2 scan of incident solar irradiance 


C4, SCAN OF TAPE8 VARIABLES 


Ifa scan of any of the TAPES variables as created by a transmittance run is required, the only 
variable treaied differently is VAR, The changes are as follows: 


IVAR 


Hi} 


O scan of HgO transmittance 
1 scan of Og transmittance 
2 scan of COs transmittance 
3 scan of CO transmittance 
4 scan of CHy transmittance 
= 5 scan of N9O transmittance 
6 scan of Og transmittance 
7 scan of NH3 transmittance 
8 scan of NO transmittance 
9 scan of NOg transmittance 
= 10 scan of SOg transmittance 
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